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GAUSSIAN ASYMPTOTICS OF DISCRETE /3-ENSEMBLES 


ALEXEI BORODIN, VADIM GORIN, AND ALICE GUIONNET 


Abstract. We introduce and study stochastic 7V-particle ensembles which are discretizations 
for general-/? log-gases of random matrix theory. The examples include random tilings, families 
of non-intersecting paths, (z, ^/-measures, etc. We prove that under technical assumptions on 
general analytic potential, the global fluctuations for such ensembles are asymptotically Gaussian 
as N —> oo. The covariance is universal and coincides with its counterpart in random matrix 
theory. 

Our main tool is an appropriate discrete version of the Schwinger-Dyson (or loop) equations, 
which originates in the work of Nekrasov and his collaborators. 
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1. INTRODUCTION 

1.1. Continuous log—gases. A general-/? log-gas is a probability distribution on /V-tuples of 
reals x\ < X 2 < ■ ■ ■ < xn with density proportional to 

N 

(1) (Xj - XiY JJexp(-W(xj)), 

1 <i<j<N i =1 

where V(x) is a continuous function called potential. For V(x) = x 2 and f3 = 1,2,4, the density 
0 describes the joint distribution of the eigenvalues of random matrices from Gaussian Or- 
thogonal/Unitary/Symplectic Ensemble; much more general potentials V(x ) are widespread and 
extensively studied in random matrix theory and bevond. see the books (Me]. |Foj. [AGZj > [02, 

m- 

Under weak assumptions on the potential, the ensembles ([Tj) exhibit a Law of Large Numbers 
as N — > oo, which means that the (random) empirical measure /xtv defined via 

1 N 

= n ^ Sxi 
%— 1 

converges (weakly, in probability) to a non-random equilibrium measure /x. For (3 = 1,2,4 and 
V(x) = x 2 , this statement dates back to the work of Wigner jWij . and /x is known in this case as 
the Wigner semicircle law. The results for generic V(x) were established much later, see [BPS] , 
|BeGu| . EH]. 

The next order asymptotics asks about global fluctuations, i.e. how the functional /xjy — /x 
behaves as N —> oo. A natural approach here is to take (sufficiently smooth) functions f(x) and 
consider the asymptotic behavior of random variables 

f(x)/j, N (dx) - J f(x)n(dx)^ , N -A oo. 

For quite general potentials V(x) the limits of 0 are Gaussian with universal covariance depend¬ 
ing only on the support of the equilibrium measure /x. In the breakthrough paper m Johansson 
proved such a statement for general /? > 0 and wide class of potentials under the assumption that 
p. has a single interval of support. Further developments have led to establishing such results for 
generic analytic potentials, see [KS] . |BoGul| . [Shdlj . [BoGu2| . Note that when the support of /x 
has several intervals one needs to be careful as an additional discrete component might appear. 
This does not happen if one deterministically fixes the filling fractions, which are the numbers of 
particles in each interval of the support. In the one-interval case the limiting covariance can be 
identified with that of a Id section of the two-dimensional Gaussian Free Field, see |B2| . |BoGo2| 
for the details. 

While for certain specific choices of potentials V(x) as well as for special values /? = 1,2,4 
there are several different methods for establishing central limit theorems for global fluctuations, 
all the developments for generic j3 and V{x) rely on the analysis of loop equations (also known 
as Schwinger-Dyson equations). These are equations for certain observables of the log-gases (flj) , 
which originated and have been widely used in physics literature, cf. [Mi] , M, El. [21 and 
references therein. More precisely, the observables of interest are the Stieltjes transforms of the 
empirical measure /x^r and the Schwinger-Dyson equations involve their cumulants. Formally 
assuming that these cumulants have a converging N^ 1 expansion as N —> oo, one can use the 
equations to derive recursively the asymptotics of the cumulants starting from the equilibrium 
measure. These recursive relations are sometimes called the topological recursion because, at least 
when p = 2, they mimic the relations between maps of different genus, see e.g. m and allow 
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Figure 1. The discrete particle systems arising in uniformly random domino 
tilings of the Aztec diamond (left panel) and uniformly random lozenge tilings 
of a hexagon (right panel). The distributions of 3 particles on both pictures have 
the form 0 with IV = 3, /? = 2 and suitable (different) potentials V{x). 

to recover the result from’t Hooft and Brezin-Itzykson-Parisi-Zuber [BlPZj showing that matrix 
integrals can be seen as generating functions of maps; for the general /? recursion see eej- 

Loop equations were introduced to the mathematical community by Johansson in m to derive 
the Gaussian behavior; his results were significantly extended in the later work |KS| . [ BoGul ]. 
|Shch| . |BoGu2| . The loop equations, or rather “their spirit”, also have further applications far 
beyond the central limit theorems for global fluctuations, e.g. they were used in the recent work 
on local universality for random matrices, see pEY| . [BKT| . We would like to emphasize an 
important distinction between the approaches of physics and mathematics literature: the latter 
operates not with formal expansions, but with converging asymptotic expansions. This requires 
a priori estimates, whose derivations rely on a set of tools different from the loop equations 
themselves. 

The Schwinger-Dyson equations for general-/? log-gases are obtained by integration by parts 
and derivation with respect to the potential. In fact, they can be derived in a much more general 
continuous setting, for instance when the underlying measure is the Haar measure on a compact Lie 
group, see e.g. [ AGZI (5.4.29)], and further used to derive the topological asymptotic expansions 
of related matrix models, see |GGM| . jGN]. For a recent application of the Schwinger-Dyson 
equations to the lattice gauge theory see [Chi] , 

l. 2. Integrable discretization of log gases. The discrete versions of the distribution |l]) (with 
Xi s living on a lattice) at /? = 2 arise in numerous problems of 2d statistical mechanics. Examples 
include random tilings (cf. Figur e [T| and ICLP], m, El, |BKMM| ), stochastic systems of non¬ 
intersecting paths (cf. KOR . IBBI )T[ ). last passage percolation (cf. |J2|). interacting particle 
systems (cf. |J3], |BF] b 

The law of large numbers for such systems (even for general values of /3) can be established by 
essentially the same methods as for the continuous ones, cf. m , du, m- However, the situation 
is drastically different for the study of the asymptotics for global fluctuations. While for some 
very specific integrable choices of the potential at /? = 2 the central limit theorem was proven (cf. 

m , m, m), until now no general approach that would work for generic /3 and V{x) existed. 
We note however the works [ HQ . Chapter 12], (DF| . |Mo| where the global Gaussian asymptotics 
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was proven for certain (very concrete) discrete general (3 probabilistic models originating in the 
representation theory. Moreover, it was not clear whether the CLT in the case of general potential 
should be the same as in the continuous case, or there should be some significant differences 
(indeed, for instance the local limits in the bulk must be different). The main technical difficulty 
lied in the absence of a nice generalization of the loop equations to the discrete setting. 

In | Ey2] , |Ey3| , Eynard proposed to interpret laws of random partitions, which include discrete 
(3 = 2 analogues of ([3]), as matrix models. One of the key points in his interpretation is to 
view discrete sums as highly oscillatory continuous integrals. Based on this identification, he 
conjectured that the asymptotic expansions in this case are described by the same Schwinger- 
Dyson equations with initial step given by the equilibrium measure of the model, see e.g |Ey2j 
Section 2.4, Section 2.7.1]. In particular, the fluctuations in the context of the central limit 
theorem should be universal. However, so far this approach has not yet progressed much beyond 
the predictions. 


The central goal of this article is to study global fluctuations of the empirical measure for 
discrete analogues of general-/? log-gases. One outcome is that indeed these fluctuations are 
universal and described by the same covariance as for their continuous counterparts. Our analysis 
is based on appropriate discrete versions of the Schwinger-Dyson equations, which (unlike in the 
continuous setting) do not appear as a direct consequence of integration by parts or perturbative 
arguments. 

The search for the discrete loop equations starts with identifying a good discrete analogue of 
the general (3 distribution 0- For that we fix a parameter 9 > 0 and a positive real-valued 
function w(x; iV)[^] Consider the probability distribution 


(3) 


Ptv(Tl, ... ,£n) 


1 ^ T(ij -£i + l)T(fy -U + 0) 

z »,ttl„ T (ej-e<Wj-e< + 1 -«) 


N 

n w&-,N) 

i— 1 


on ordered iV-tuples i\ < £2 < ■ ■ ■ < £n, such that £$ = A i + 9i and Ai < A 2 ■■■ < \n are integers. 
We refer to £fi s as the positions of N particles. Let us note that if 9 1, then £i do not sit on 

the fixed lattice. 

Note that if 9 = 1, then ([3]) has the same form as ([Tj) with (3 = 2. Similarly, 9 = 1/2 leads to ([I]) 
with (3 = l|^J More generally, if we set £i = Lxi, then as L — > 00 , the ratio of Gamma-functions in 
(|3j) behaves as {£j — £i) 2e and mimics 0 with (3 = 29. 

The most important reason to view (Kh) as a correct integrable discretization of the continuous 
log-gas 0 is the following observation which is the starting point of the results in the present 
paper. 


Theorem 1.1 (Nekrasov’s equation). Consider the probability distribution ([3]), and assume that 

w(x] N) = <j>tr( x ) 
w(x — 1; N) </>-( x ) 

and define 


(4) fliv(0 = <MO-E P, 


N 


n 1 


_ 2 — 1 


Z-ti 


+ </>+(0-Ep, 


N 


_ 2 — 1 


n>+ 


e - u -1 


If are holomorphic in a domain JCl N c C, then so is Rjy(Z)- Moreover, if are 

polynomials of degree at most d, then so is Rn(Z)- 


^w{x\ N) should decay at least as |x| (2JV 2)0 1 const with const > 0 as |®| —¥ 00 . 
2 However, observe that the lattice where particles sit at 9 = 1/2 is not Z. 
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Theorem |1.1| is essentially due to Nekrasov and his collaborators, as it is a variant of similar 
statements in [N] , m, |NSj . For the proof see Theorem |4.1| below. 


It is reasonable to ask how one could guess the form of the measure ([3]), which would lead to 
Theorem 0 The integrability properties of such measures, including the product of the ratios 
of Gamma functions as in ([3]), can be traced to their connections to representation theory and 
symmetric functions. The same products of Gamma functions appear in the evaluation formulas 
for Jack symmetric polynomials (cf. |Mal Chapter VI, Section 10]) and in problems of asymptotic 
representation theory (cf. |Q2j h Another trace of integrability is the existence of discrete Selberg 
integrals, which are evaluation formulas for the partition function Zj\t in ([3]) for special choices of 
the weight w(x-,N). Examples of such evaluations can be found in [ 02t Section 2], [ GSl Section 
2.2], The latter reference also explains the degeneration to the conventional Selberg integral, 
which is the computation of the normalization constants for continuous log-gases ([!]) with specific 
choices of the potential V(x), see [Me, Chapter 17], [El Chapter 4], On the other hand, for a 
“naive” discrete version, when one takes the same formula ([!]) as a definition of a discrete log-gas, 
we are not aware of the existence of similar evaluation formulas (outside (3 = 2, when (jTJ) and Q 
coincide). 


1.3. Global fluctuations of discrete log gases: one cut case. We proceed to our results on 
global fluctuations for the distributions ^ as N —> oo. For that we need to postulate how the 
weight w(x;N) changes as N -A oo. Our methods work for a variety of possibilities here, but to 
simplify the exposition we stick to the following assumption in this section: 

(5) w(x-N) =exp (~ NV (j^)) , 

where V{z) is an analytic function of real argument z such that for |z| large enough, V(z) is 
monotone and satisfies 


(6) V(z) > c ln(| 2 :|), where c > 26. 

Our first result is the law of large numbers for the empirical measures hn defined via 


N 


hN = 


J^J2 6 £i/N, 

i =1 


(£±,... ,£ n ) is IPw — distributed. 


Theorem 1.2. There exists a deterministic absolutely continuous compactly supported probability 
measure p(x)dx with 0 < /j,(x ) < 6~ x , such that /u/v converges to p,{x)dx as N -A oo, in the sense 
that for any compactly supported Lipshitz function f(x) the following convergence in probability 
holds: 

lim / f(x)p, N {dx) = / f(x)n(x)dx. 

N ^°° Jr Jr 

In fact, we prove a more general statement where, in particular, V(x) does not have to be 
analytic, see Theorems 5.3 and 10.1 below. The measure p,(x)dx is the equilibrium measure, and 
it can be found as a solution to a variational problem. At 6 = 1, Theorem |5.3| reduces to results of 
.13j. [J4j, |Fej . For general values of 6, additional arguments are required, and we present them. 

Note the condition /r(x) < 9~ 1 , which is not present in the continuous log-gases, and arises 
from the fact that the minimal distance between adjacent particles in ([3]) is at least 6. This is a 
specific feature of the discrete models. 

At this moment we need to make certain assumptions on the equilibrium measure. A frane^of 
n{x) is a maximal interval (a,/3) such that 0 < n(x) < on ( a, (3 ). From the random matrix 


3 We follow the terminology from BKMM i. 
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literature (cf. |BDE] . [Shch] . [BoGu2j ) one expects that the global fluctuations are qualitatively 
different depending on whether p(x) has one or more bands. Here we stick to the one band case. 

Theorem 1.3. Assume ( [5|, (JFi] ), that p,(x) has a unique band (a_,a+), and that (technical) 
Assumption^ from Section \3.I\ holds. Take any m> 1 bounded analytic functions /i ,..., f m on 
M. Then the m random variables 

N 

C fj = J2(m) - N fj(ti))> C^i) ■■■An) is P n -distributed, 

i=l 


converge (in distribution and in the sense of moments) to centered Gaussian random variables 


with explicit covariance depending on a± and 6, and given in (127) below. In particular, if 
fi(x) = (z,i — x ) _1 , Zi G C \ R, i = 1,..., m, then 


(7) Jim 


1 ~ 


ZiZj — 2 (a_ + a+)(z+Zj) + a+a_ 
yj( z i - a-)(zi - a + )\/( z j ~ a~)(zj - a+) 


Remark 1.4. We prove below in Theorem |10.1 that with exponentially high probability all particles 
ti are inside an interval [—DN,DN]. Further, it is enough to assume in Theorem 1.3 that /,; are 
analytic only in [-D, D], Moreover, we believe (but do not prove) that the analyticity assumption 
can be replaced by sufficient smoothness. 

Remark 1.5. The second order asymptotic expansion of the mean Ep^ J)(fi)^J can be also 

analyzed, cf. Theorem 6.1 with k = 1, m = 0. 


The proof of Theorem |1.3| is a combination of Theorem 7.1 and Theorem 10.1 in the main text. 

The technical Assumption [4] from Section 3.2 is a statement that a certain function produced 
from the equilibrium measure n(x)dx has no zeros. In Section 9.3 we show that this assumption 
always holds when V (x) is a convex function. A somewhat similar assumption appears in the work 
of Johansson m on the central limit theorem for global fluctuations of continuous log-gases. 


Let us emphasize that the limiting covariance in Theorem 1.3 depends only on the support 
of the equilibrium measure, but it is not sensitive to other features. The same phenomenon is 
known in the random matrix setting. Moreover, comparing ([T]) with expressions in |J11 Theorem 
4.2], [PEi Chapter 3], we conclude that the covariance is precisely the same as for the continuous 
log-gases. Thus, the discreteness of the model is invisible on the level of the central limit theorem. 

The covariance in Theorem P can be related to that of a section of the 2d Gaussian Free 
Field in the upper half-plane with Dirichlet boundary conditions. One way to predict that is by 
noticing that sections in lozenge tilings models with several specific boundary conditions yield 
distributions of the form ([3]) with 6 = 1, see Section 9.2 for one example. On the other hand, 
there exist several results on the appearance of the Gaussian Free Field in the asymptotics for 
global fluctuations in lozenge tilings, cf. E, ML n |BuGo2]. 

Subsequent work |KS| . |BoGul) for continuous log-gases showed that Johansson’s technical as¬ 
sumption on the equilibrium measure can be relaxed and replaced by certain weaker assumptions, 
which hold for generic analytic potentials. We hope that a similar thing can also be done in the 
present discrete setting, but this would require further investigations. 


1.4. Weight supported on several finite intervals. In several applications the weight w(x; N ) 
is not defined on the whole real line, but instead it is supported by a union of several disjoint 
intervals ( ai(N),bi(N )), i = In other words, the particles ti are now confined to the 

union of these intervals (J^i (oi(N),bi(N)). For instance, this happens in tilings models, see 
Section 19.21 
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At this point we have an additional choice: One could either fix the filling fractions rii(N), 
which are the numbers of particles ii in each of the intervals {afiN), bfiN)), or not. In the present 
paper we stick to the former and fix m(N). 

We refer to Sections [3] and [5] below for details of the exact assumptions that we impose on 
the weight w(x ; N ), intervals (ai(N), bi(N)), and filling fractions rii(N). Here we will only briefly 
summarize the results obtained in such a framework. 


The first result is the law of large numbers, which is the exact analogue of Theorem 1.2 


see 


Theorem 5.3 for the details. As before, for the asymptotics of the global fluctuations we need 


to have some information about the bands of the equilibrium measure p,{x)dx. We assume that 
there is one band per interval. 


Theorem 1.6 (Theorem 7.1). Assume that all the data specifying the model satisfies Assumptions 
011 in Section [iO| below. In particular, the equilibrium measure p,(x)dx has k bands (aj,/3j), 
i = 1,..., k, one per interval of the support of the model. Take m > 1 functions fi(z), ..., f m (z), 
which are analytic in a neighborhood of U^ =1 [a/(A r ) /N,b{(N)/N] for large N. Then as N —y oo 
the joint moments of the m random variables 

,... ,£n) is Pn- distributed, 
approximate those of centered Gaussian random variables with asymptotic covariance depending 


N 

c-u =EWF*)-Ep„/j(C). 

i= 1 


only on 9, a ±,..., a*,, fii, ■ ■ ■, fik and given by (116). 


For k = 1, the limiting covariance coincides with that of Theorem 1.3 For k > 1 it is no longer 


expressible through elementary functions, e.g. for k = 2 elliptic functions appear in the formulas. 
We again emphasize that the covariance depend only on the end-points of the bands and is not 
sensitive to other features of the equilibrium measure. Furthermore, the covariance is the same 
as for the continuous log-gases, cf. fShch] , |BoGu2| . 

Similarly to the k = 1 case, one could try to identify the covariance with that of a canonically 
defined random field. We believe that there should be a link to sections of the Gaussian Free 
Field in a domain with k — 1 holes, but we postpone this discussion till a future publication. 

1.5. Organization of the paper. In Section [2] we present our approach and results: Nekrasov’s 
equation, Law of Large Numbers and Cent ral Limit Theorem — in the simplest yet non-trivial 
case when 9 = 1 and 




in Theorem 


1.1 


are linear functions. The measure P/v in this case 


is known as the Krawtchouk orthogonal polynomial ensemble; it appeared and was studied in 
numerous previous articles. In particular, the Central Limit Theorem for global fluctuations in 
this specific case can be also obtained by other methods: see mn], \m, and |BuCo2( 5 

In Sections [3]-[8] we explain a much more general framework, in which the same ideas work and 
lead to the Law of Large Numbers and Central Limit Theorem. In these sections 9 is an arbitrary 
positive number, and the weight w(x\ N ) and functions are general (subject to certain technical 
assumptions). 

In Section [9] we specialize the general framework to certain specific examples, which include 
lozenge tilings and (z, ic)-measures from asymptotic representation theory. We explain how all 
the technical assumptions are checked in all these examples. 


Finally, Section 10 explains how the case of the infinite support of w(x]N ) can be reduced 
to the case of the bounded support (which is studied in Sections [3]-[8]) through large deviations 
estimates. 


4 We remark that while the approaches of these articles admit various generalizations, to the best of our knowledge, 
they do not extend further in the direction of the present paper. 
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2. Toy EXAMPLE: 6 = 1 BINOMIAL WEIGHT 

The aim of this section is to describe our method in the simplest, yet non-trivial case of the 
binomial weight and 6 = 1. The resulting IV-particle ensemble is known as the Krawtchouk 
orthogonal polynomial ensemble. It appeared in the literature before in the connection with 
uniformly random domino tilings of the Aztec diamond (cf. m ), last passage percolation (cf. 
[J2]), stochastic systems of non-intersecting paths (cf. [KOR] - (BBDTj), and with representation 
theory of the infinite-dimensional unitary group U( oo) (cf. [ BQ2 . Section 5], [B1L Section 4]). 

In the subsequent sections we will extend our approach to more general cases, but the method¬ 
ology and many key ideas remain the same. 


Fix two integers 0 < N < M and consider the space Wat of A-tuples of integers 

0 < i\ < £2 < • • • < &n — M. 

We define a probability distribution Pjy on Wtv through 


( 8 ) 


Pat(^i, ..., In) = 


1 


N 


Z(N, M) 


n (<i-4) 2 n 


M 


l<i<j<7V i— 1 

We remark that the partition function Z(N, M ) is explicitly known in this case: 

IV -1 


Z(N, M) = 2 iv(m-tv+ i )( M! )IV j-j 


J 


n 


j =0 


(M-j)f 


However, in a generic situation there are typically no simple formulas for the partition function 
and we are not going to use its explicit form. 

Our analysis of the distribution Fn is based on the following 6=1 version of Theorem 4.1 This 
is essentially due to [Nl . m, [NS] and we will call the main statement the Nekrasov’s equation. 
It should be viewed as a discrete space analogue of the Schwinger-Dyson (also known as “loop”) 
equations. 

Proposition 2.1. Let P^v be a probability distribution on N-tuples 0 < l\ < ■ ■ ■ < £n ^ M such 
that 

1 N 

F N (£ 1 ,...,£ N ) = - ;Q {l, - W (4), 

7=1 


where 


l<i<j<N 

w(x) _ 4 >n( x ) 


^(a?- 1 ) <M®)’ 

and are analytic functions in a complex neighborhood of [— 1, M + 1] which are positive on 

[0, M] and satisfy cj)fr(M + 1) = = 0. Define 


(9) Rn(0 = ■ IEpj> 


N 


n- 


_i=l 




+ <&(£)■ Ep, 


N 


.7=1 


n 1+ 


t-b-i 
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Then Rn{Q is analytic in the same complex neighborhood of[—l,M+l\. If(f N (£) are polynomials 
of degree at most d, then Rn {£) is also a polynomial of degree at most d. 

Proof. The possible singularities of Rn(£) are simple poles at points m = 0,1, 2,, M. Let us 
compute a residue at such a point. 

The expectation Ep^ in © is a sum over all possible configurations (JL\,,£n) £ Wat. Such a 
configuration contributes to the residue at m if either ^ = raor^ + l = m for some i = 1,..., N. 

We separately analyze the contributions appearing from each i = 1 ,,N, which we now fix. 
Given a particle configuration £ = (fi, •.., £n), let i + denote the configuration with ith coordinate 
increased by 1 and let l~ denote the configuration with ith coordinate decreased by 1. Note that, 
in principle l + (similarly £~) might fail to be in Wtv, as the coordinates might coincide. However, 
in this case the formula for Pjv(-f + ) still applies and gives zerc0 

The contribution to the residue of Rn at m. E {0 ,,M}, arising from the ith coordinate of 
(h ,... ,l N ) is 


( 10 ) 


(t> N (rn)P N (£i,.. 

^GWjvl t-i=m 






m — In 


+ 

feWjv|L=m-l 


m — £j — 1 


But using the definition of Pjv we see that 


F n (£i, .. .£i-i,m,£ i+ i,... ,£ N ) JJ M 


1 


m — in 


= (t>^(m)F N (£ 1 ,. • .ii-i,m - l,i i+ i, • ■. ,^jv) JJ f 1 + m _ ^ _ A 

j+i V 3 ' 


We conclude that for each pL, the terms with i = pi and i = pi + (or £ = pi~ and £ = pi) in the first 
and second sum in (10) cancel out and the total residue is zero. 

For the polynomiality statement it suffices to notice that if 0^(£) are polynomials of degree at 
most d, then Rn(£,) is an entire function which grows as 0(£ rf ) as £ —> oo. Hence, by Liouville’s 
theorem Rn(£) is a polynomial. □ 


Note that for the distribution ([8]) the functions 4>%{x) can be chosen to be linear and we set 

A-t \ x a+( \ M + 1 x 

( n ) M x ) = n’ <w) = —^ 

Thus, in a sense, ([8]) is one of the simplest possible distribution in the framework of Proposition 

E31 

In this section we aim to study the asymptotics of the distributions Pjy as IV —> oo. The 
parameter M will also depend on N. We fix m > 1 and set M = [nrN"J. 


2.1. Law of Large Numbers. As N —> oo, a certain law of large numbers holds for the measures 
P N- Let us introduce a random probability measure pijy on M via 


( 12 ) 


1 


N 


^ n = n^2 6 

i =1 


(£i, ... ,£n) is P/v-distributed. 


J This is where we need the condition of + 1) = <j> N ( 0) = 0, which translates into w(M + 1) = w(— 1) = 0. 
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The measure /iy is often referred to as the empirical measure of the point configuration £ 1 ,..., £y. 
Note that our definitions imply the condition £i + \ — ii > 1, which shows that for any interval 
[p, q], its measure is bounded from above by (q — p + 1 /N). 

The idea for the proof of the law of large numbers for /ry as N —» oo is to establish the large 
deviations principle for the measures Py, which would show that the measure is concentrated on 
£i which maximize the probability density ([8]). This was done rigorously in m, [J4] , see 
also Section 2.4 below for some details. The explicit formula for limy^oo/ry is then obtained as 
a solution of a variational problem, this solution for our weight was found earlier in jDSl Example 
4.2], The developments of these articles are summarized in the following proposition. 

Proposition 2.2. The measures pn converge (weakly, in probability) to a deterministic absolutely 
continuous measure p m (x)dx, which is called equilibrium measure. Form > 2, the density of the 
measure p m (x)dx is 


hm(x) = < 

and for 1 < m < 2, 


1 

— arccot 
7r 


m — 2 


Lo, 


/hnO) = < 


1 

— arccot 

7T 


2\/m — 1 — (x — m/2 y 


m — 2 


2\/m — 1 — (x — m/2)^ 


1, 

0, 


\x — t)| < V'R - l, 


\x — > \/>Tl - 1, 


\x — | < ynx - 1, 

^ > \x — ?r| > -\/m — 1, 

I _ m I s, m 
rn 9. — •?. 1 


where arccot(y) is the inverse cotangent function. 

A convenient way of working with the equilibrium measure /i m (x) is through its Stieltjes trans¬ 
form G m (z) defined through 

1 p m (x)dx 


(13) 


G m {z) = 


z — X 


Observe that (74) makes sense for all z outside the support of p. m (x), and G m (z ) is holomorphic 
there. Further, as z —> oo, we have G^z) ~ z~*. An explicit formula for G m (z ) can be readily 


extracted from Proposition 2.3 below. 


We introduce two functions 


R m {z) = zey:p(-G m (z)) + (m - z) exp(G m (z)), 

Qm{z) = zexp(-G m (z)) - (m - z) exp(G m (^)). 
for z£C\[0,m] and consider their analytic continuations. 


Proposition 2.3. For any m > 1, 

R m (z) = m - 2, Q m (z) = 2 y/(z - m/2) 2 - (m - 1), 
with the branch of the square root chosen so that Q m (z) ~ 2 z as z —> oo. 


is, 


Proof. We first observe that R m (z 
limit of Rjy(Nz) in Proposition 2. 
implies lirn^oo R m (z) = m — 2, and thus R m (z) = m — 2. 
Q m and R m implies 


in fact, a linear polynomial. Indeed, 


this is the N oo 

,-i 


for the measure Py. The fact that G m (z) ~ z " as z —> oo, 
Further, observe that the definition of 


Qm(z) 2 - R m (z) 2 = -4z(m - z). 
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Therefore, 


Qm(z) 2 = 4 z 2 — 4zm + (tn — 2) 2 = 4(z — m/2) 2 — 4(m — 1). □ 


2 . 2 . Second order expansion. Define the Stieltjes transform Gjy(z) of the prelimit empirical 
measure ( 12 ) through 

/ °c 1 1 N ^ 

^ _ m(dx) = — ^2 _ n / t, r i (h, ■■■An) is P/v-distributed. 

-oo Z ^ Z tj/iV 

We aim to study how Gn(z) approximate G m (z) as N — > oo. For that we introduce a deformed 
version of the same function. Take 2k parameters t = (t\,.. . , t/-), v = (yi,...,Vk) such that 
v a + t a — y 7 ^ 0 for all a = 1 , ,k and all y E [0,m] and let the deformed distribution P^ v be 
defined through 


(15) P % v (h,...,£ N ) = 


1 


N 


Z(N,M- t,v) 


n »- i i > 2 n 




Z— 1 _ 


If k = 0, then P^ v = P/v is the undeformed measure. In general, P)y V may be a complex-valued 
measure. The normalizing constant Z(N, M;t,v) in (15) is chosen so that the total mass of P^) v 
is 1, i.e. X^eWjv CO = L Let us note that the numbers t a are always chosen small enough, 
which guarantees that Z(N, M;t, v ) ^ 0. 

Observe that P^r V satisfies the assumptions of Proposition | 

k 


M 

£i 


7)t. V 


a= 1 


np+ 


tn 


Va ~ ii/N 


2.1 


with 


(16) 


(17) 


0 iU®) = 


M + 1 x 
N 


N 

k 


a =1 


nb-p/ /---I • 


N N 


a =1 


As above, we set M = [mAj, and omit it from the notations. We also define 


Clearly, 


^N (0 = ( l mN \ + 1 - mN ) + Y 

a =1 

</>n( N z ) = 


m — z 


v n - z 


0 n( z ) = I] 


d— 1 


t a + V a - Z 


+ 0 | jv^ 


n (w - z) ( t a +v a ~z), 


a =1 


<f> N (Nz) = (z + 


0 n( z ) 

A 


+ 0 l iV^ 


n ( U “ “ ~) ( ta + v a ~ z) . 


a= 1 


We further define /ijy V as the empirical distribution of P^ v and set 


t,v 


/ CX) 1 

- Vn ( dx) 

-oo Z-X 

Define 



2=1 


1 

z-ti/N’ 


(£i,... ,£n) is P^ v -distributed. 


(19) A G n (z) = N(G n (z) - G m (z )). 

Note that we often omit the dependence on t, v from the notations, to keep them concise. 

The definition of the deformed measure P*i/ V is motivated by the following observation. It was 
used before in the related context in random matrix theory, cf. [IMi] , m - 
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Lemma 2.4. For any k > 1, the kth mixed derivative 

< 20) atl9t2 ... att E -y( AG ”W) 

£ a = 0 , l<a<fc 

is the jomf cumulant of k + 1 random variables NGn(u), NGn(vi), ..., NGj\r(vk ) with respect 
to the measure Pat. 

Proof. We first note that the joint cumulants are invariant under addition of constants. Therefore, 
we can replace NGn(u ) by A Gn(u) for the purpose of computing the cumulants. 

Further recall that one way to define the joint cumulant of k + 1 bounded random variables 
Co, • • • ,6 is through 

hi f E exp f ^2 U£i 


^fc+i 


dtodti... dtk 


o=o 


to=ti=-=tk=0 


Taking the derivative with respect to to explicitly, we can rewrite this also as 

Qk E (Co exp (j2i=i ) 


dti ... dt k 


E ( ex P (Ei=l ti& 


N 


ti=-=t k =0 

Setting Co = NAGn{u), Ci = NGn( vf ), i = 1 ,,k and observing that 

P + -‘ 

i—1 

we get the desired statement. 


ex.p(tNG N (v)) = JJ ( 1 + 


v - li/N 


+ 0(t 2 ), t —> 0, 


□ 


Theorem 2.5. Fix k = 0,1, ... and choose k complex numbers v = (ui, ... ,Vk) E (C \ [0,m]) fc . 
Then, as N —> oo, 


( 21 ) 


d k 


dt\ ■ ■ ■ dtk R N 
dz 


E p t,v (A Gjv(h)) 


d k 

= o(l) + - 

t o =0, l<a<fc dt\---dtk 


1 


27ri ■ 2^(w — m/2) 2 — (m — 1) 


5 


7[0,i 


(u-z) 


^-( z )e- G ^+^(z)e G ^ + V(z - m/2) 2 - (m - l)-^G m (z) 


_ t a =0, l<a<fc 


where 7[o,m] a simple positively-oriented contour enclosing the segment [0,m] (the points u and 
v\,... ,Vk are outside the contour). The remainder o(l) is uniform over u,v i,.. •,Vk in compact 
subsets of the unbounded component of C \ 7[o,m] • The k = 0 case is that we take no derivatives 
in pll). 


Remark 2.6. The only cases, where the right-hand side of (21) is meaningful are k = 0,1, since 
its N —> oo limit is zero for k > 2. Indeed, only depends on t a , and moreover it is a sum 
of functions of single variables t a , a = 1, therefore, all mixed partial derivatives vanish. 

However, we present Theorem 2.5 in this way, since that’s the form which appears in our proofs. 
For k = 1 we will compute the limit of (21) below in Theorem 2.8 For k = 0 and generic m the 
right-hand side of (21) has no limit, as QmAj + 1 — miV) in the definition of oscillates. 


Remark 2.7. It might seem a bit unexpected that for k = 0 the right-hand side of (21) does not 
vanish when 6 = 1. Indeed, in the context of random matrices this corresponds to fd = 2 case, 
where the mean is known to vanish (cf. |Jlj ). For our model the non-zero mean can be traced back 
to two features. First, due to discreteness we can not adjust M = [mAj so that ([mA] +1 — mA) 
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oo 


becomes zero for all N. Second, the logarithm of the weight, In has a non-trivial N 

asymptotic expansion, which affects the result. 

Before providing a proof of Theorem 2.5 let us give its important corollaries. Note that if we 
set k = 0 in Theorem 2.5, then it gives the limit behavior for Ep^ AGn(z). In fact, it also gives 
the Central Limit Theorem for Gjy(z). 

Theorem 2.8. The random field N^Gn(z) — z G C \ [0,m] , converges as N —> oo 

(in the sense of joint moments, uniformly in z in compact subsets of C \ [0,m]/ to a centered 
complex Gaussian random field with second moment 

(22) lim N 2 (~E Pn (G n (u)G n (v)) -~E Pn G n (u)E Pn G n (v)) =: C{u,v) 

uv — ^(a_ + a+)(u + v) + a + a_ \ 


1 


1 - 


2(u - v ) 2 y y/{u — a-)(u - a+)y/(v - aJ){v - a+) J ’ 

where a± = ™ =t Vm — 1. 

Remark 2.9. Note that the covariance C(u,v) has no singularity at u = v, since the right-hand 
side of (22) has a finite u —> v limit. 


Remark 2.10. Since Gn(u) = Gn(u), the formula (22) is sufficient for determining the asymptotic 
covariance of the random field Gn{u)- 


Proof of Theorem \2.8[ We start from the result of Lemma 2.4 for the joint cumulant of k + 1 
random variables NGn(u), NGn(v\), ..., NG^ivk)- In particular, if k = 1, then we get the 
covariance of NGn(u ) and NGn(v i). 

For k > 1, differentiating © we see that the result vanishes as IV — > oo, see Remark |2.6[ This 
implies the asymptotic Gaussianity of the random field N^Gisr(z) — Ep jv Gat(z)J. 

In the case k = 1, differentiating (21) we see that the covariance given by (20) is (recall that u 
and v\ lie outside the integration contour) 


(23) 


1 


2-7ri • 2 \J(u — m/2) 2 — (m — 1) 


dz 


/ 7[0,n 


(u — z) \ (ui — z) 2 


0 G m (;z) 


dz 


2vri • 4^(11-m/2) 2 - (m - 1) J 7[0 , (z - «)(«i ~ ^) 2 


Rm(z) + Qm(z ) 


The term with R m (z) integrates to 0, as it has no singularities inside 7[o. m ]- The term with Q m is 
computed as the sum of the residues at z = u and at z = v \, which gives the desired covariance 
formula. □ 


Theorem 2.8 implies the central limit theorem for general analytic linear statistics. 

Corollary 2.11. Let the distribution P n be given by ©> with M = [mlVj . Take k real valued 
functions f\ (z),, fk{z) on [0,m], which can be extended to holomorphic functions in a complex 
neighborhood B of [0,m]. Then as N —> oo the k random variables 

N 

C fj = J2(m) - e Pjv/j(^))> {h, ---An) is Pa redistributed, 

i =1 

converge in the sense of moments to centered Gaussian random variables with covariance 


(24) 


lim Ep Ct.Cf. = t 

A^oo rN L h ( 2 7ri) 2 


fi(u)fj(v)C(u,v)dudv, 


,m] , m' 
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where 7[o m ] is a positively oriented contour in B which encloses [0,m], and C(u,v) is given by 

(22b- 


Remark 2.12. The covariance (24) has the same form as for random matrices and log-gases in the 
one cut regime. It depends only on the restrictions of functions fj onto the interval [a_, a+] and 
can be rewritten in several other equivalent forms, cf. m Theorem 4.2], m Chapter 3], [ACZ , 
Section 4.3.3]. 


Proof of Corollary \2.11 


Observe that 


£/ = 


27ri hw 


f(z ) {G n (z) - Ep N G N (z)) dz. 


Therefore, all the moments of Cf are obtained from the centered moments of Gn(z) by integration. 
Since the latter converge uniformly in z on the integration contour, so do the former. It remains 
to use the fact that the integrals of jointly Gaussian random variables are also Gaussian. □ 


Remark 2.13. Similarly to Corollary 2.11, Theorem |2.5| can be used to obtain the first two terms 
in the IV —» oo asymptotic expansion of Ep^ YliLx /(4) for functions / holomorphic in a neigh¬ 
borhood of [0, m]. 


The rest of this section is devoted to the proof of Theorem 2.5 


2.3. Heuristic argument for Theorem 2.5 , In this section we present a sketch of the proof for 


Theorem 2.5 in which we omit crucial bounds on remainders in the asymptotic formulas. These 
bounds will be established further on. 


We start from the statement of Proposition 2.1 for the measures P^- v . Making the change of 
variables £ = Nz, we have for z outside [0,m] 


N 


(25) n 1 


i= 1 


£-4 


= exp 


N ( 1 


s. 2=1 


N z-tifN 


- exp ( -C»(z) + — ■ qzG n (z) + O ( A ,. 

= exp ( -<?„(*) - iAC„( 3 ) + ^ (e) + O (T 


where the remainder is uniform over z in compact subsets of C \ [0, m]. Similarly, 


N 


( 26 ) n !+ 


2=1 


Z-£i- i 


= exp I In I 1 + — • 


N z-li/N-l/N 


= exp ( G m (e) + fAGrj(z) - T , JLg n ( z) + O (T 
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Recalling the definition of R m (z), we conclude that the function Rjy(Nz ) from Proposition 2.1 
can be written in the following form 


(27) Rn{Nz) = JJ(u a - z)(t a + v a - z) 


a =1 


Rm{z ) 


+ ze Gm(2:) E p t,v ^exp ^J^AGn(z)\ - 1^ + (in - z)e Gm E p t,v ^ 


exp ( - — AG N (z) ) - 1 


, -G m (z) , G»(z) 


+ 


ze Gm ( 2 ) fid 


N \2dz 


n^Z G m(z) + 


(m — z)e GmG ' > / 19 


N 


2 dz Gm{z ' 


+ ° { N 


where appeared from the second terms in 1/N expansion of (fr^ from ©, ©• We fur¬ 
ther want to simplify the expression in the second line of (27), by replacing e h — 1 by h under 
expectations. For that we note a simple inequality, which we will use with n = 2: 

71—1 


(28) 




b? 




e - y — 


j =0 


<\h\ n e W , he C, n = 1, 2,... 


We will later establish in Section 2.5 that uniformly in z in compact subsets of C \ [0, m] we have 


(29) 


E. 


K 


G n (z) 


exp 


j^AGn(z) 


= °( jy ) 1 N ^°°- 


We therefore can rewrite (27) as 
k 

(30) R n (Nz) = Y]_(v a -z)(t a + v a -z) ■ 


a =1 


R m (z) + (A Gn(z)) 


+ 


M z ) c ~G m (z) . 'I’nW c G m (z) , ^ Gm{z) ( 1 9 G , A (M^ (_\d_ 

N e + N 6 + N \2dz Um[Z) ) + N V 2 8z 


'dz Gm ^ z) 




where the estimate of the remainder is uniform over 2 in compact subsets of C \ [0, m]. 
Let us now fix u outside the contour 7[o, m ] enclosing the interval [0, m], divide ( f30j ) by 

k 

2 ■ 2tt\ ■ (u — z) ■ JJ(u a - z)(t a + V a - z) 


a=1 


and integrate over 7[o, m ]- Since both Rn(Nz) and R m (z) are holomorphic inside the contour, the 
integrals of the corresponding terms vanish. From the rest we get, with the help of Proposition 

[231 


< 31 > h 


V( z ~ a -)( z ~ a +) 


7[0,n 


U — Z 

1 

27ri 


Ept.v (AGn(z)) dz 

r N 


7[0,. 


dz 

2 (u — z) 


ipfr(z)e~ Gm{z) + ^(z)e GmG) 


+ ze~ G ^ Q^G m (z)) +(ra-z)e G ^ (-^G m (z))) +o(l). 
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For the first line of (31), note that E p t,v (AGjy(z)) is analytic outside the contour of integration 

and decays as 1/z 2 when z —>■ oo. Therefore, we can compute the integral as (minus) the residue 
at z = u, which is 

y/(u — a-)(u - a+) • E p t,v (A Gjy(u)) . 


This leads, via differentiation in t’s, to the desired formula (21) of Theorem 2.5 


Let us point out the parts of the above argument that are not yet rigorous and whose justifi¬ 


cation is necessary to complete the proof of Theorem 2.5 


We need to show that the bound (29) is valid. 

We need to prove that all the remainders remain small when we differentiate with respect 
to variables t a , a = 1,..., k. 


2.4. A weak a priori estimate. The following lemma is a first step for establishing the desired 
estimates. 


Lemma 2.14. Fix a positive integer n and take a compact set A C C \ [0,m]. Then for every 
e > 0 there exists a constant C > 0 such that for every zi ,..., z n € A and every N = 1,2,... we 
have 


(32) 


Ei 


'Pjv 




2—1 


< c ■ N n ^ +£ \ 


with AGn(z) defined by (19) with k = 0. 

Remark 2.15. We will show in the next section that C ■ N n ^ + A in the righthand side of (32) can 
be replaced by a constant. However, to produce such a sharp estimate we need to start from this 
weaker one. 


The proof of Lemma 2.14 that we now present is similar to the argument of |BoOu2l Lemma 
3.4.5], which follows the ideas of jMMSj . 

Before we proceed we need to recall the characterization of the equilibrium measure y m . Intro¬ 
duce the functional I[p\ of a measure /ion [0, m] via 


(33) 


If A = JJ ln \x - y\dfj,(x)dia(y) - j V (x)dfi( a 


0<T,y<m 

xAy 


where 

V(x) = x ln(x) + (m — x) ln(m — x). 

Note that if the measure p has atoms, then it is important to exclude the diagonal in the integral 
( |33] ), i.e. integrate only over x / y. 

The variational characterization of the measure / u m (x)da0 (see [33], m, ra yields that 
p m (x)dx is the unique minimizer of I[fi\ among the probability measures of density at most 

1. 

Define the functions F m (x) through 


r m 

(34) F m (x) = 2 / ln \x — t\fj, m (t)dt — V(x). 

Jo 

Then varying the functional at /j m , one proves (cf. |DSj . [|Fe], [ST]) that there exists a real 
number / such that F m (x) — / = 0onS’ = {0<x<m|0< /i m (x) < 1}, F m (x) — / < 0 on the 
complement of the support of ji m . and F m (x) — / > 0 when the density /x m (®) is equal to 1. 


^Throughout the paper the density of a measure y is denoted y(x). 
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Now take any two compactly supported absolutely continuous probability measures with uni¬ 
formly bounded densities v(dx) = u{x)dx and p(dx) = p(x)dx and define V(y, p) through 

(35) V 2 (v,p) = -[ [ ln|x - y\{u{x) - p{x)){v{y) - p(y))dxdy. 

J M J M 

There is an alternative formula for T>(v(x), p(x)) in terms of Fourier transforms, cf. [BeGuj : 

2 


(36) 


V(y,p) = ^ J^e itx {u(x) - p(x))da 


dt. 


Fix a parameter p > 2 and let pn denote the convolution of the empirical measure p^ given 
by (12) with uniform measure on the interval [0, N~ p ], 

Proposition 2.16. There exists CgK such that for all 7 > 0 and all N > 1 we have 

F N (v(p N ,Pm) > 7 ) < exp(CiVln(iV) - 7 2 iV 2 ). 

Proof. Observe that for every iV-tuple 0 < l\ < I 2 < ■ ■ ■ < In < M we have, using Stirling’s 
formula for factorials, 


(37) P N (h,...,£ N ) = 

where 


exp^2iV(lV- l)ln(lV) + N 2 l[mes[£i,... ,i N ]] +0(NI 

Z(N, |mlVj) 


N 


mes[£i, ■ ■ -Jn] = 


i —1 


N 


Let us obtain a lower bound for the partition function Z(N, LmiVj) in (37). For that let Xp 
i = 1,..., N be quantiles of p m defined through 

r x i i — 1/2 

J p m {x)dx = ———, 1 < i < N. 

Since p m (x) < 1, Xi + 1 — Xi > l/N and therefore, the numbers , 1 < i < N are all distinct. 

We can then write 

(38) Z{N, |mlVj) >exp(2N(N- l)ln(JV) + JV 2 /[mes[LlVxiJ,..., [Nx n \]] +0(N In 
We claim that ( |38[ ) can be transformed into 

exp(2N(N - l)ln(N) + N 2 1 [p m ] + 0(Nln(N))^J . 

Indeed, for the double-integral part of /[•] we write using the monotonicity of logarithm 


X> 

i<j 


\Nxj\ I NxA \ v ( 1 

Xj-X, + - 


N 


N 


i<j 


< N ' 


^ r x j +1 r x i 
, , ■ Jxj JXi -1 


i<j *3 


In ( t - s + — ) p m (t)pm(s ) dt ds + 0(N In(IV)) 


= N 2 


ln(t — s)p m (t)p m (s) dt ds + 0(N In(IV)), 


S<t 


and similarly for the opposite inequality. For the single-integral part of /[•] we have 
[Nxi 


nE v 


1=1 


N 


JV 2 ^ 


1,+I V(t) + (t- FTa 1 v"(«(t)) Wt)* + 0 (N ), 


— 9 J Xi 


i =2 


N 
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where n(t) is a point inside [xi,xjv]- Observe that |W(k(£))| = 0(ln(lV)), and thus 

N , , Ar , s N—2 


nJ2 v 


i =1 


L Nxi 

N 


n 2 E 


Ei + l 


i=2 


V ( t ) + ^x i+ i - Xi + ) 0(ln(iV)) ) p m (t)dt + O(iV) 

N-2 


= N 2 J\(t)^ m (t)dt+N 2 ° {ll ^ N)) J2 (*i+i - = N 2 J\(t)y m (t)dt+0(Nln(N)). 

The next step is to replace mes[fd,... ,£n] i n (37) by its convolution with the uniform measure 
on [0, N~ p ], that we denote riles[0 ,... ,£jy], For that take two independent random variables u, 
u uniformly distributed on [0, N~ p ], where p > 2 as above. Then 


(39) l[mes[£i,... ,£ N ]] = E Ui £ 


/ / In \x — y + u — {t|mes[^i,... ,£iv](dx)mes[^i,... ,0v](dy) 

/ o Jo 

rm 

— ~E U V(x + u)m.es[£i,... ,£]y\(dx) 

Jo 


1 

= /[mes[£i,... ,^jv]J + / In |it - u|mes[7i,... ,£x](dx) 

Jo 


+ 


In 


x^y 


1 + 


N 
u — u 


x-y 


mes[7i,...,£jv](d®)mes[£i,... ,£]y]{dy) 


+ E U J (V(x + u) - V(x))mes[^i,... ,£ N ](dx) = J[mes[^i,... ,£ N ]] + O • 

We conclude that there exists a constant C such that 

¥ n (£i, ■ ■ ■ ,£n) < exp(CiVhi(iV))exp^iV 2 (/[mes[7i,...,£Ar]] -/[/x m ]) 

Further, completing the square we get 
(40) l[mes[£i ,.. .,£ N ]] ~ 7[/hn] 

pm 

=-'D 2 (p m ,rnes[£ 1 ,...,£ N ]) + / F m (x)(mes[£i,... ,£ N ] -/j, m ){dx) 

Jo 

pm 

= -V 2 mes[£i ,... ,£ N ]) + / (F m (x) - f)(mes[£i,... ,£ N ] - p m )(dx), 


with F m {x) and / defined in (J34J) and directly below that formula, respectively. Let us analyze 
the last term in (40). On the set S = {x | 0 < /i m (x) < 1} the function F m {x) — / vanishes. On 
the complement of the support of /r m (x) we have F m (x) — / < 0 and the corresponding part of the 
last integral in (40) is negative. Finally, on the set S' = {x | p m {x) = 1 }, we have F m (x) — / > 0. 
Since all the points £i/N are at least 1/N apart, for any a such that [a, a + 1 /N) C S' we have 


(mes[0, ...,£n\~ y m ){[a,a + 1 /N}) < 0. 


Thus, 


( 41 ) 


ra+l/N 


< 


(F m (x) - /)(mes[4, ...,£ N ]~ p m )(dx) 

f-a+l/N 2 

(F m (x) - F m (a))(rnes[0, ■ • • ,£n] ~ dm){dx) < — sup \F m (x) - F m (y)\ 


N 


x,y£S’ 
\x-y\<l/N 
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Observe that the last sup is at most 0(ln(N)/N). Now partition S' into segments of the form 
[a, a + 1/1V] and note that for boundary segments (which have to be shorter than 1/N) the bound 
for the integral of the form (41) is still valid as F m (x) is equal to / in one of the end-points of 
such segment. Summing the bounds over all segments, we get a bound on the integral over S'. It 
follows that as N —> oo 

J (F m (x) - f)(mes[£i,...,£ N ] - y m )(dx) < O . 

Therefore, we finally obtain 

Fjv(4, ■■■,£n)< exp(C' N \n(N)) exp( -N 2 V 2 (mes[h,.. .,In],L m)Y 


Since the total number of IV-tuples 0 < i\ < .. .In < is (^at + 1 ) = exp(0(iVki(lV))), the 

proof is complete. □ 

Corollary 2.17. For a compactly supported Lipschitz function g define 


11/2 


/ OO 

e isx g(x)da 

-OO 


2 \ >/2 
ds 


I Lip = sup 

x^y 


g(x) - g(y) 


x-y 


Fix any p > 2. Then there exists C £ M such that for all 7 > 0, all N >1 and all g we have 
(42) F n 


pm pm 

/ g(x)/j, N (dx) - / g{x)n m (dx) 

/ 0 Jo 


> 7119111 / 2 +^) 


< exp ( CNln(N) — 


7 2 Ar 2 


Proof. We have 

i pm 

g(x)fi N (dx) - / g(x)y m (dx 

Jo 

pm 

/ g{x)fi N {d. 

Jo 


< 


x — 


g(x)fi N (dx) 


+ 


i pm 

g(x)fi N (dx) - / g(x)y m (da 

Jo 


The first term is bounded by ^ p ip and corresponds to such term in (42). Therefore, it remains 
to work out the second term. Since scalar products are preserved under Fourier transform, with 
the notation 


/ OO 

e lsx (f(x)da 

-OO 


we write the Plancherel formula (note that g(x) and (/2 at(x) — fi m (x)) are bounded and belong to 
Z / 1 [0, m] n L 2 [0,m]) and use the Cauchy-Schwartz inequality 


g(x)(fl N (x) - y m (x))dx 


f 


AaV) -Am(9 

vW 


dt 


< llslli/21/ r lw(t) = v^llslli / 2 ®(w, 9 m). 


It remains to use Proposition |2. 16 


□ 


Proof of Lemma 2. If. Choose a small g > 0 and take an infinitely differentiable function h(x), 
whose support is inside [— 17 , tn + 7 ] and such that h(x) = 1 for 0 < x < m. 
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Since both fi^ and // m are supported on [0,m], we can replace l/(z — x) in the definition of 
Gn(z) and G m (z) by a nice smooth compactly supported function h( x)/(z — x) without changing 
Gn(z) and G m (z). Now choose 7 = q ■ N ~ l / 2+£ , q > 0, in Corollary 2.17 and note that for this 

00 . This readily implies the 

□ 


choice the right-hand side of (42) still exponentially decays as N 
bound of Lemma 12.141 


2.5. Self-improving estimates and the proof of Theorem 2.5| . The aim of this section is 
to finish the proof of Theorem 2.5. The key part is establishing the following statement. 


Proposition 2.18. For any k > 1 and any u±,... ,Uk £ C \ [0, m], the joint moments 


Ep^ XI I AGjv(«a)| 


a =1 


are uniformly (in N and in u\,... ,Uk in compact subsets of C \ [0, mJJ bounded. 


The idea for getting such estimates is to start from Lemma |2.14| and then recursively feed the 
existing estimates into the argument of Section |2.3| to obtain the stronger ones. This is very 
similar to the argument of jBoGuh, Section 4.3]. 


Proof of Proposition 2.18. Fix n = 0,1,... and v = (i?i,..., v n ). Arguing as in Lemma 2.4 
prove that for any bounded random variable £ we have 


we 


(43) 


dP 


dt\ • • • dt r 


E, 


a) 


= M c (Z,NG n (v 1 ),...,NG n (v„)), 


t\= — =t n = 0 

where M c is the joint cumulant. Since the cumulants are unchanged under shifts, we can also 


replace NG]y(v a ) by AG]sr(v a ) in (43). 


We aim to differentiate the formulas of Section |2.3| with respect to t a at t a = 0. For that we 
need to examine the remainders. The remainder o(lV _1 ) in ( |27| ) comes from three sources: from 
the expansions (25), (26); from the IV -1 expansion of (jr^(z)', from the replacement of -^Gn(z) 
by -^G m (z). The remainder can then be written as a sum corresponding to these three sources: 
(44)“ 

n 

Y[(v a -z)(t a + v a -z) 

a=1 


j^K p t^^ N (z) + j^c(z;t, v)E p t^v£)y(2:) + (^,n{ z )—AGn(z) 


where £jv(*) are random variables, which are bounded uniformly in N, in z belonging 

to compact subsets of C \ [0,m], and which do not depend on t, v. The function c(z;t,v) arises 
from the large N expansions of ©, ©, and it is uniformly bounded in z belonging to compact 
subsets of C \ {ui,..., u n ; ii + ui;... t n + v n }. The dependence on z is holomorphic in all the 
terms. 

Further, when we pass from (27) to (30), we need to expand Kp N (exp (-^AG/v^)) — l). Since 
jfAG]\r(z) is a bounded random variable, we can use the Taylor expansion to get 


(45) E„t;v 


exp 


N AG N (z)j lj E p v 


— AG N (z)^j + 2 E p^ v 
1 


A G N (z] 


+ A 

6 r N 


A G n (z] 


+ ... 


Thus, after reconstructing the remainders, (31) is replaced by 
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(46) 

1 

47ri 


R m (z) 


Qm{z ), 




u — z 


• Q m (z)E p t,v (AG«(z))+^E F t,v (AG/v(^) 2 ) + ^ Z E p t,v (AG/v(*) 3 )+. • • 


+ 


N 


Ept;v^jv(^) + c(^;;t, v)E p t;v^(^) + E p t;v ( (z)—AG n (z) )) 


I dz 


1 

27ri 


7[0,t 


dz 

2 (u — z) 


■ (^-(z)e- G ^ +^+(z)e G ^ + ^ ) ^G m (z) 


We now differentiate (46b with respect to all t a at t a = 0. For all the terms involving ran¬ 
dom variables we use (43) to rewrite the result as a cumulant. Note that when we differentiate 
c(z; t, v)E]pt;v£jy( 2 :), we need to apply the Leibnitz rule and therefore to differentiate each of the 

two factors and get a sum. We also compute the integral of ;A^Q m (z)E p t,v (AGn(z)) as minus 
the residue at z = u to get 

Qm{u) 


(47) 


2 -M c (AG/v(u), ACtjv(^i)) • • • j AG]\r(v n )') + 

jf ^ (^M c ((AG ff (.)) 2 , AGjv^i), • • •, A G N (v n )) 


+ Q^M c ((AG n (z)) 3 , A G n ( Vi ), ..., A G N (v n )) + ... j dz 


—^—;M c U n (z) + £x(z)-^-AG n (z), AG n (vi), AG N (v n ))dz 
J'Vm m i U Z V C>Z ' 


1 1 

N 4ni ,/ 7 , u — z 


1 1 
N 47ri 


E 


d r 


dt\ ■ ■ ■ dt n 


1 

4-7ri 


Ac{i,...,n} J7 [°. m ] 
dz 

7[0,m] (« “ Z ) 


u — z 


-c A (z-,t,v)M c ((;' N (z),AG N (v a ),a <E A)dz 


■( ^^)e“ Gm(z) +^^)e Gin W + ^^-^G m (^) 


_ t a = 0,1 <a<n 


where the fifth line is the result of the differentiation of c(z; t, v)E p t ; v£j( r (,z), i.e. CA(z;t,v) is the 
mixed derivative of c (z\ t, v) with respect to t a , a G {1,..., n} \ A. Let us analyze the resulting 
expression (47) as N — > oo. At this moment all we need from the right-hand side of (47) is that 
it is 0(1) as N —> oo. 

Consider the infinite sum over growing powers of AGn(z) in (47). All the terms starting from 
the H th one can be combined into 


(48) 


Me ( E (Q/ R )h ( ^[ Z } r AG/vM, • ■ •, AG N (v n ) ) , 

\h=H ' / 


where ( Q/R)h is either Q m (z) or R m (z) depending on the parity of h. Expanding the cumulants 
in terms of centered moments and using Holder’s inequality, we observe that if £ is a bounded 
random variable, then for any fixed k > 1 there is a constant C k such that 

k 

(49) \M c ((,AG N (vi),...,AG N (v k ))\ < C k [J^IAO/v^) " E[AG/v(^)]| 2fc ] ^ 

i =1 
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Observe that the random variables Gi\r(v a ) are uniformly bounded, same is true about ( Q/R)h 
(for z and v a in compact subsets of C \ [0, m]), and 


c = © ( Q/R ), 


h=H 


(A G n (z)Y 
2h\N h ~ l 


satisfies (uniformly over z in compact subsets of C \ [0, m]) the point-wise bound 


|C| < const • N 


A G n (z) 


N 


H 


exp 


AG'tv(z) 


N 


Thus, we can bound (48) by 
(50) 


const -N n+1 sup 


N 


IEp A 


AGjv(z) 


N 


2 H 


exp 


2AG n (z) 


N 


< const-./V n+1 sup 


N 


E P A 


A G n (z) 


N 


2 if N 


where in the last inequality we used the uniform boundness of AG ^( U ) , anc [ SU p^ ^- 0 taken 
over any set which includes both 7[o, m ] and all points v a , a = 1,n. 

We can use Lemma 2.14 to bound the expectation in (50) and conclude that if H > 3 n, then 


(48) is o(l). Therefore, renaming u into vo, (47) is rewritten as 
Qm(llo) 


(51) 


-M c (AGtv(l , o), AGat(ui), ..., AGt\t( fn))T 


+ 


1 


1 


4 - 7 ri L, , vq — z l A-/ jyh l 

' 7[0,m] u \h =2 


3 n \ 

E • ^M c ((AGv(z)) k , a G N ( V1 ), A G N (v n )) dz 


+ 


- — l —M c U N (z) + ^(z)~AG N {z), AG n ( Vi ), AG N {v n j)dz 

J 7rn V 0 Z ' GZ 


1 1 

N 47ri J-noM v ° - z 

1 1 

+ -•-71 E 


N 47ri *—' , vn — z 

Ac{l,...,n} 7l °’ m 


ca(z; t, v)M c (C(v(z), AG N (v a ), a e A)dz = 0(1), 


where ( Q/R)h is either Q m (z) or R m (z) depending on the parity of h, and we replaced u by vq. 
At this moment we claim that © for all n = 0,1, 2 ... together with the bound of Lemma 


2.14 implies Proposition 2.18 Indeed, take any two disjoint compact sets W,VcC\[0,m], which 


are invariant under conjugation, and suppose that 7[o,ml C U . Expanding cumulants in terms of 
centered moments and using Lemma 2.14 we obtain 


1 1 


(52) sup . 

Vo,—,V„&V A 47T1 


t -^—M c U N {z)+^(z)^-AG N (z),AG N (y 1 ),...,AG N (v n ))dz 

i /v rn ml u 0 - z V dz ) 


= O 


n+ l)(l/2+e)-l 


(53) sup — • — / - 

vo,~-,Vn£V A 47T1 A J^ n} T 1[0 m] V 0 - Z 


C a(z; t, v)M c (g N (z), AG N (v a ), a <E A)dz 

= O ^Ar n l 1 / 2+e ) _1 j 
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(54) 


sup 

voy,v n ev 


4-7ri 


T[0,m] 


' ^!^M c ({AG N {z)) h , AGiv(ui),..., A G N (v n ))^j dz 


= O 


n) (1 /2-\-e) + 1 —h 


One might be cautious about (52), as it involves the derivative of AGn(z) (instead of AGn(z) 
itself as in Lemma 2.14), yet for analytic functions a uniform bound for a function implies a 
uniform bound for its derivative, therefore, the bound is valid. 


Qm(^o) 


We plug (|52|), (|53|), (|54|) into (|5l|) to get 
(55) 


sup 

vo,...,v n eV 


M c [AGn(vo), AGat(ui), ..., A Gn(vti)) 


+ 


3 n 




lh=2 


ri)(l/2-\-e)+l—h 


+ O ^lv(" +1 )( 1 / 2+E )- 1 ) + o =0(1), 


which implies that for each n = 0,1,2... 

(56) sup | M c (AG n (vq), ..., AO/vK-i))! = 0(N^ + i)(i/2+e)-(i/ 2 - e )) + 0 ( 1 ), 

vo,vi,...,v n eV 

where the remainder depends on the choice of the compact set V C C \ [0, m]. As centered 
moments are linear combinations of products of joint cumulants, we deduce from (56) that for all 
A: = 1,2,..., 


(57) 


sup E 

vi,...,v k eV 


Q(AOjv(uo) - E[AGiv(u 0 )]) 


.a =1 


= o( N W*+*)-W-*)) + 0(1). 


Combining n = 0 version of (56) with (57) we finally conclude that 

k 


(58) 


sup E 

v\,...,Vk£V 


11 A G N (v a 


a =1 


= o 


+0(1), A: = 1,2,. 


If k is even, then we can choose v\ = ■ ■ ■ = v k / 2 = v and v k / 2 +i = ■ ■ ■ = v k = v in (58) and get 


(59) 


supE 

vev 


|AO/vM|* 


= O 


e)-(l/2-e)\ +0(1), As = 2,4,6,_ 


The estimate for odd k is reduced to the even k case with the use of Jensen’s inequality in the 
form 


E(£ fc ) < (E (t k+1 )) 

which leads to the following bound for all k: 


k 

k +1 




(60) 


supE 

vev 


\AG N (v)\ l 


= O 


, k > 0, 

e) —^ t( 1/2 - £ ))+0(1), As = 1,2,3,.... 


We finally use the Holder inequality to get 
k 


(61) sup E 

vi,...,v k eV 


11 \AG N (v a )\ 


.a =1 


< supE 
vev 


\AG N (v)\ k = O (V fc ( 1/2+£ )-fcT t(i/ 2 - £ )^ + o(l) 

2+£ )-i/ 6)+0(i), As = 1,2,3,... ; 


where we silently assumed that e > 0 is small enough for the last equality to hold. 
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At this moment we can iterate the argument. Expanding cumulants in terms of moments, we 


deduce from (61) the following three bounds: 


(62) sup ' 7 -t —-— M c (f >N (z)+ f , f r (z)-^-AG N (z),AG N (v 1 ),...,AG N (v n ))dz 

vo,—,v n eV N 47T1 J 7 j 0 m ] Vq z \ oz J 

= o ^Ar( n + 1 )( 1 / 2+£ )- 1 - 1 / 6 ^ , 


(63) sup ^ / —-— 

„ 0r ..,, ne v N 47T1 ( ^ / 7[0 , m] V 0 - z 


c A (z-,t,v)M c (g N (z), AG N (v a ),a <E A)dz 

= O (V n(1 / 2+£) “ 1-1 / 6 ) 


(64) 


sup 

vo,...,v n eV 


4-7ri 


T[0,m] 


vq — z \N h 


r • { ~m^M c ((AG N (z)) h , A G N ( Vl ),A G N (v n ))\ dz 


= O 


(^N {h+ 


n)(l/2+e)+l— h— 1/6 


The only difference between (52), (53), (54) and (62), (63), ( fcij) is t hat the degree of N in the 
bound decreased by 1/6 (which is because the bound of Lemma 2.14 is replaced by (|61[)). 

We then plug (l62j), ([631) , (l64| into (l5lj) to get 


(65) sup 9^.M c (AG N (y Q ), AG N ( Vl ),..., AG N (v n )) 

vo,—,v n eV 2 


+ 


3 n 

^ () (/H-n)(l/2+s)+l-fc-l/6^ _j q ^jy(ra+l)(l/2+e)-l-l/6\ _|_ Q — 0(1). 


lh=2 


In the same way as (55) implied (61), the bound (65) implies 

k 


(66) sup E 

v 1 ,...,v k eV 


n |AGjvK)| 


a= 1 


= o 


2+ £ )-4t(1/2- £ )-1/6^ +0 ( 1 ) 


= o 


-2' 1 / 6 )+0( 1), k = 1,2,3,.... 


Repeating the same argument m — 2 more times, we improve (66) to 

k 


(67) sup E 

vi,...,v k eV 


II IAGatK) 


a= 1 


= o 


e) - ro/6 )+0(l), k = 1,2,3, 


Since m is arbitrary, this implies 


( 68 ) 


sup E 

vi,...,v k eA 


n \*G N (v a )\ 


.a= 1 


= 0(1), A: = 1,2, 


and finishes the proof of Proposition 2.18 


□ 


Proof of Theorem 2.5 Take (47) and observe that the bounds of Section 2.5 imply that all the 
terms except for the second line and the last line are negligible as N —> oo. The second line of 
(47) is precisely the left-hand side of (21), while the last line of (47) is the right-hand side of 

(2l). " □ 









































GAUSSIAN ASYMPTOTICS OF DISCRETE /3-ENSEMBLES 


25 


3. General setup 


3.1. Definition of the system. Our next goal is to generalize the arguments of the previous 
section to a much more general setting of a multi-cut fixed filling fractions model with fixed 
parameter 6 > 0 and general weight w(x). Informally, we want to consider measures of the form 

T{ij - e, + i)r(< 3 - - e, + o) N 


n 




r (£j - -£i + 1 - 6 ) 


n w^n) 


i=i 


on ordered IV-tuples l\ < £2 < ■ ■ ■ < £n referred to as positions of N particles and satisfying 
two additional constraints. First, the particles are separated into k groups, and particles in each 
group must belong to its own interval of the real line. Second, if ith and ( i + l)st particles are in 
the same group, then £i + 1 — li G {6, 9 + 1, 9 + 2, ... }. 

For instance, if we have a single group (k = 1), then after defining A/ through li = \i + 6i , the 
constraint boils to down to saying that all A/ are integers and they satisfy Ai > A2 > • • • > Ajv- 


Let us give precise definitions for general k. The model depends on an integer parameter 
N = 1,2,... and amounts to fixing for each N a probability distribution on certain IV-point 
subsets of M. 

We fix an integer k = 1,2,..., whose meaning is the number of segments in the support of the 
measure. For each JV = 1,2,... we take k integers m(N ),..., rik(N), such that Yli =1 n *(-^0 = N 
and k disjoint intervals (ai(N),b\(N)), ..., ( a,k(N),bk(N )) of the real line ordered from left to 
right Q 

We assume that bi(N) + 6 < a,i + i(N) for i = 1,..., k — 1. The numbers cii(N ), bi(N) must also 
satisfy the conditions 

(69) bi(N) - Orii(N) - a*(IV) G Z. 

Further, the number n/(fV) counts the number of the particles in the ith interval; to make this 
statement precise we define the sets of indices Ij C {1,... ,1V}, j = 1 ,,k, via 

f ' 7_1 J 

Ij = li€z\ ^2 n m (N) <i n m(N) 

\ m=l m =1 

We also set if and I~ to be the maximal and minimal elements of L, respectively. 

3 3 J 



Definition 3.1. The state space W 0 N consists of IV-tuples £1 < £2 < ■ ■ ■ < £n such that for each 

j = 1, ■ • • ,k: 

(1) If i = IT, then li — a,i(N ) G Z>o- 

(2) If i = If, then bi(N) - l t e Z >0 . 

(3) If i € Ij, but i f Ij~, then £i + 1 — £i € {0, 6 + 1, 6 + 2,... }. 


Note that the conditions of Definition 3.1 imply that for every i from Ij, we have li G [a*(IV) + 
1, bi(N) — 1]. An example of a configuration from W 0 N is shown in Figure [2J 

We also take a positive weight function w(x;N) for x G ujL 1 [a/(A r ) + 1 ,bi(N) — 1] and define 
a probability measure Pat on W e N given by 


(70) 


Piv(*i, • • • ,In) 


1 n riij-ii + i^iij-ii + d) 

Z N ^v^ij-iiWj-ii + i-e) 


N 

l\w(li-,N), 

2=1 


where Zn is a normalizing constant which we will refer to as the partition function. 


7 n 

For a generalization to the case of infinite support see Sections 


9.3 


9.4 
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int 9 + int int 


int 


int 


<2.1 (3) (-i 


M3) 


02 ( 3 ) 


6 2 (3) 


m(3)=2 particles 


77 , 2 ( 3 )=! particle 


Figure 2. The state space for N = 3, k = 2. Numbers int indicate various 
(possibly different) nonnegative integers. 


3.2. List of regularity assumptions. Our ultimate goal is to study the asymptotics of the 
measures Pat as N — > 00 . For that we need to assume that the weights w(x;N), as well as all 
other data specified in Section 3T depend on N in a regular way. Let us present all the technical 
assumptions that we impose on the data. 


Assumption 1. We require that for each i = 1, ..., k, as N —> 00 

ai(N) = Ndi + 0(ln(lV)), bi(N) = Nbt + 0(ln(A)), 
di < 61 < 02 < • • • < dk < bk- 

We require that w(x] N) in the intervals [a* (IV) + 1, bi(N) — 1], i = 1,..., k, has the form 

w(x;N) = exp (~NV n (^)) 

for a function Vjy that is continuous in the intervals [a* (N) + 1, bi(N ) — 1] and such that 
(71) Vh(u) = V(u) + o(^P-') 

uniformly over x = u ■ N in the intervals [oj(iV) +1, bi(N ) — 1]. The function V ( u ) is differentiable 
and the following bound holds for a constant C > 0 

k 

1 + y~](|ln(tt- dj)\ + I ln(u - Sj)|) . 

1=1 

Remark 3.2. We believe that the assumption on the remainders can be weakened with minor 
changes in all the further statements and proofs. However, we do not pursue this direction due 
to lack of natural examples. 

For the filling fractions nt(N) we make a weaker assumption, as it might be important for 
future applications, cf. [BoGu2j. 

Assumption 2. There exists a constant C > 0 such that for all 1 < i < k and all large enough 
N we have 

c < r h^l <e -i { j ) ._d i )-c. 

Note that our definition of the state space W e N implies that < 6~ 1 (bi — d{) + o(l). 

Introduce the notation 

Ui 

— = n i; i = l 

Note that the numbers hi still depend on N. However, we will hide this dependence from our 
notations. It is important here that all the limiting values as well as all the remainders in what 
follows will be uniform over hi satisfying Assumption [2] 


(72) 


\V'(u)\ < C 


The next two assumptions deal with analytic properties of the weight w{x\ N ) and the equilib¬ 
rium measure p. We fix an open set M.n C C, such that uf =1 [oi(IV), bi(N)] C Mn- 
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Assumption 3. There exist a pair of analytic in x £ At at functions (f>~^(x), <t>^{x) such that 

w{x\N) = 4>n{x) 
w{x -1 ;N) </>Jj(xY 

Moreover, 

^ (v) + (v) + ° (vy) 

uniformly over x/N in compact subsets of an open set A4 ac which contains the union of the 
intervals [a/, bi], i = 1 All the aforementioned functions are holomorphic in At at and 

functions ip ^ are uniformly bounded as N —> oo. 

Remark 3.3. In the case when Vn(x) in Assumption [l] is smooth and uniformly converges to V 
together with its derivative in a neighborhood of a point x, we have 

<73) exp (-= Pw>- 

Indeed, in this case, 


4> + (x) 

4>-(x) 


lim ^(-NA’N) = lim exp(-N(VN{x) — Vn{x — 1/N)) = 

3V->oo w(Nx — 1; N) oo v v ’ v 111 


= lim exp —N 

TV-s-oo \ 


/ Kv(y)dy ) ) = exp 
Jx-l/N J J 



Recall that the equilibrium measure p with the density p(x) encodes the Law of Large Numbers 
for F n stated in Theorem |1.2| and discussed in more details in Section [5} A convenient way of 
working with the equilibrium measure is through its Stieltjes transform G^z) defined through 


(74) 



p(x)dx 
z — x 


The following two functions R^z), Q^(z) are important for our asymptotic study, cf. Section 

EZQ 

(75) := exp (-dG^z)) + 4> + (z) exp^G^z)), 

(76) Qfi(z) ■■= 4>~(z) exp (-OG^z)) - 4> + {z) exp (OG^z)). 

We explain in Section [5] that R^(z) is analytic, while Q^z) is a branch of a two-valued analytic 
function which is the square-root of a function holomorphic in -Mat. An important technical 
ingredient of our method is a restriction on its zeros, as is summarized in the following assumption. 


Assumption 4. We require that for each large enough N and corresponding Q^ (which de¬ 
pends on N through the filling fractions hi in the definition of p), there exists a function H(z) 
holomorphic in At at and numbers {a/,/3i}( ; L 1 such that 

• hi < a.i < ^ < bi, i = 1 ,..., k; 

• Qu( z ) = H ( z ) n \/{ z - a i)( z - a), 

i =1 

where the branch of the square root is such that ^{z — a.i)(z — /?/) ~ z when z oo; 

• H(z) ^ 0 for all z 6 U/=i[dj, &/]• 
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We remark that Assumption [4] does not describe a generic case. In particular, it implies that 
there is precisely one interval of support of n{x) in each interval [di,bi\. The authors are not 
aware of simple ways to check such property by examining the potential V(x). Nevertheless, 
many natural models arising in the applications satisfy Assumption [4j We demonstrate this 
general principle by considering several examples in Section [9j 

Finally, we need a simple vanishing assumption. It is convenient to work with it, yet we later 
show in Section [8] how it can be relaxed. 

Assumption 5. For all * = 1,, k, we have </>iy((a,;(A r ) + 1) = ^(bfiN)) = 0. 


4. Nekrasov’s equation 

The main tool for our study of the probability distributions Pat from the last section is a 
statement, which is essentially due to Nekrasov fNj, m, fNSj- Its affine and ^-versions are given 
in the following two theorems. 


Theorem 4.1. Let P^r be a distribution on N-tuples {£\, . . . , In) E 
Suppose that 

(77) w(x;N ) = fi + N {x) 

w(x-l;N) <fQ(x)’ 

and for all i = 1,..., k, we have <fifij(ai(N) + 1) = <fi^(bi(N)) = 0. Define 

9 


/ 9 n as in the previous section. 


(78) i2;v(0 = <M0-E p. 


' N 

n 

,i =1 


i - 


S-ti 


+ </>+(£)• E P , 


N 


. 1=1 


n i+ 


9 


Z-ti-1 


If 0^(0 are holomorphic in a domain M.n C C, then so is Rn(€)- Moreover, if are 

polynomials of degree at most d, then so is i?jv(£). 

Proof. The possible singularities of Rn{£,) are simple poles arising from the denominator of the 
expression under expectation Ep^ in (78). Let us compute a residue at such a pole m. 

The expectation Ep^ in (78) is a sum over all (£\,... ,£n) £ W e N . Such a configuration con¬ 
tributes to the residue if £i = rn or = rn + 1 for some i = 1 ,... ,N. 

We separately analyze the contributions appearing from each i = 1,..., N, which we now fix. 
According to definitions, the possible values for are {A, A + l, A + 2 ,..., B} for certain A and B. 
Given a particle configuration £ = {£\,..., £n), let £ + denote the configuration with Ah coordinate 
increased by 1, and let £~ denote the configuration with Ah coordinate decreased by 1. Note that 
in principle, £ + (similarly £~) may fail to be in W e N . However, in this case the formula for Pjv(^ + ) 
still applies, but gives zero. 

Let us explain how the weight (70) of a configuration (G, • • • ,-^v) changes when one coordinate 
is changed from = x to £i = x — 1, i.e. we compute the ratio of the weights at £i = x — 1 and 
at £i = x (all other coordinates are unchanged). The double product over i < j in (70) produces 
factors (we denote £j = r here) 

T(r — x + l)T(r — x + 9) T(r — x + l)T(r — x + 2 — 9) (r — x)(r — x + 1 — 9) 


[> ' :h r(r - x)r(r - x + 1 - 9) 

if i < j, and the factor 

r(x — r + l)T(x — r + 9) 


(80) 


T(r — x + 2)T(r — x + 1 + 8) (r — x + l)(r — x + 9)' 
r(x — r — l)r(x — r — 9) (x — r)(x — r + 9 — 1) 


r(x — r)r(x — r + 1 — 9) T(x — r)T(x — r + 9 — 1) (x — r — l)(x — r — 9) ’ 


if i > j. Note that (79) and (80) are two forms of the same rational expression. 
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Now take m G {A, A + 1,..., B}. The contribution to the residue of Rn at z = m, arising from 
the ith coordinate of configurations (£],... ,£n) is 


(81) -6 Y <f N (m)F N (£ 1 ,...J N )]j(l--^] 

■feW® | ti=m jAi ^ ^ 


+ 6 


Y <t>N(m)F N (ei, ■ ■ ■ ,MlI ( 1 + 

j¥=i 


i—m—1 


m — £j — 1 


Note the difference between two summation sets in (81). If £i is not the smallest particle in an 
interval [dh(N) + l,bh(N) — 1], then the first sum contains the terms with \ = m — 0, while 
the second one does not. However, each such term in the first sum is actually zero. Also if £/ is 
the smallest particle and m = A = a,h(N ) + 1, then the second sum is empty as £i is never m. 
But we also know that (p7(rn) = 0^(a/ l (A r ) + 1) = 0 and the first sum vanishes as well. Similar 
considerations apply to two cases whether £i is the largest particle or not. We conclude, that it 
suffices to study the case when there is one-to-one correspondence between terms of two sums in 

481b- 


Using (jT9j) , (|80j) , and ([77]) we see that 
<j>~ N {vn) Pjv(£i, •. .£i-i,m,£i+i,... ,£ N ) ^1 


m — £a 


= <t>N(m)F N (£i,.. .£i-i,m - l,£ i+ i, ...,£n) U + m _ ^ _ 1 j 

j^i ^ j / 


We conclude that for each p, the terms with £ = p and £ = p + (or £ = p~ and £ = p) in the first 
and second sum in (81) cancel out and the total residue is zero. 

For the polynomiality statement it suffices to notice that if are polynomials of degree at 

most d, then Rn(0 is an entire function which grows as 0{f d ) as £ —> oo. Hence, by Liouville’s 
theorem Rn(£) is a polynomial. □ 


The proof of Theorem |4.1| reveals that it admits a natural (/-deformation. Recall the definition 
of (/-Gamma function T g : 

r g (x) = (i - qf~ x {q ' q)co 


(q x -,q )c 


where 


(a;q) oc = f|(l 


aq 


n =0 


In the same framework of Section [3j define a (/-deformation of Pat through 


(82) P» r (( 1 ,...,< w ) = d r 


n x 

N l<i<j<N 


-e^-U) T q {tj-£i + l)T q {£i-li + °) 


N 




Theorem 4.2. Let P^ be a distribution on N-tuples (£\ ,... ,£n) G W 0 n as above. Suppose that 

w(x;N) (/>n{x) 


w(x — 1; N) (t> N (xy 

and for all i = 1,..., k, we have 4>Jj(ai(N) + 1) = = 0. Define 


(83) i4(0 = <M0-E p® 


N 


n * 


_i=1 


1 - q^~ e 
1 — q£~ £i 


+ ^n(0 ' 


N 


n« 


-i=l 


1 _ q i-L-i+e 


si — q 
2 


1 - q ^- 1 
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If r/»±(0 are holomorphic in a domain M. N c C, then so is R q N (£). Moreover, if are 

polynomials of degree at most d, then so is Rn(£). 

Proof. The proof is the same as for Theorem |4.1| The only new ingredient is a g-deformation of 
(79), (80), which now reads 


(84) q 


-o(r-x) r g( r - x + 1 ) r g( r -x + 9) e(r _ x+1) Tq(r - a + 1 )T q (r - x + 2 - 9) 

T,(r - x)T q [r - x+l-9) T q (r - x + 2)T q (r - x + 1 + 9) 


= q 


,(1 ~q r ~ x )(l -q r ~ x + 1 - e ) 

— - qr-x+O'j 


(i — q 


and 


(85) q 


-O(x-r) T g( x - r + ^(g -r + 9) g( x _ r _i ) rqr(x -r- l)r q (x -r-0) 


T g (x — r)T q (x — r + 1 — 9) 


T g (x — r)T q (x — r + 9 — 1) 

9 (l-q x - r )(l-q x - r+e - 1 ) 

(1 - q x ~ r ~ l )(l - q x -r - e )' 


= q 


One readily recognizes the same expression in (84) and (85). We further conclude that 

el — q m ~ e i- e 


(f N (m)P q N (£i,.. .£i-i,m,£ i+ u ... ,i N ) [q* — 




qin—tj 


f el — a m ~ 

= <^(m)P^i,...^_i,m- l,£ i+1 ,...,£ N ) JJ (q~2 - _ m _ £j _ 1 J 


and therefore there is a cancelation of the poles. 

5. Law of Large Numbers 


□ 


5.1. Limit shape. Assumptions [l] and [2] guarantee that a certain law of large numbers for the 
measures Pjy holds as N —> oo. To state it we introduce a random probability measure pn on M 
via 

i N 

( 86 ) 


1 


i— 1 


Pn = ^2 5 ( ~k ) ’ (^i, ■ • •, In) is Pjv-distributed. 


The measure pn is often referred to as the empirical measure of point configuration £\,... ,£jy, 
cf. (12). Note that our definitions imply the condition £{ + i — £i > 9, which shows that for any 
interval \p, q], its pn -measure is bounded from above by 9~ 1 (q — p + iV _1 ). 

The first order asymptotic behavior of measures pjy can be understood through a variational 
problem. For a probability measure p supported on Uj[dj,6j], define 


(87) 


I v [p\= 9 I In \x - y\p(dx)p(dy) - 
Ix^y 


V (x)p(dx). 


Lemma 5.1. Let 0 be the set of absolutely continuous probability measures p(x)dx supported on 
U i[di,bi\, whose density is between 0 and 9~ x and such that 

Pi 

/ p{x)dx = hi, 1 < i < k, 

J CLi 

where 0 < hi < 9~ l (bi — af), i = 1,..., k, are such that = 1- Then the functional Iy has 

a unique maximum p{x)dx on 0. 


Remark 5.2. The maximizer of Iy is called the equilibrium measure. 
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Proof of Lemma 5.1 . Let us equip © with weak topology, i.e. the topology of pointwise conver¬ 
gence for the distribution functions. Then 0 is compact. 

Observe that the functional ly is continuous in the weak topology on 0; here it is crucial that 
the measures in 0 have density between 0 and Q^ 1 . Therefore, ly attains its maximum on 0. It 
remains to show that such a maximum is unique. 

For that we note that ly is strictly concave, i.e. for any p, p' G 0 and any 0 < t < 1 

I v [tp + (1 - t)p] > tly [p\ + (1 - t)l v [p]. 

Indeed, the linear part of ly is concave by the definition. The quadratic part is negatively-definite 


due to the formula (36) for it, therefore, it is strictly concave. 


Since a strictly concave functional cannot have more than one maximum on a convex set, we 


are done. 


□ 


Theorem 5.3. Under Assumptions^ an <*0 the random measures pjy concentrate (in probability) 
near p(x)dx of Lemma \ 5.1\ More precisely, for each Lipshitz function f{x) defined in a real 
neighborhood o/ujL-Jd/, bj\, and each e > 0 the random variables 


iv 1 / 2 - 


f 


f(x)p N (dx) - / f(x)p(x)dx 


converge to 0 in probability and in the sense of moments. 

Remark 5.4. Since p(x) depends on hi, and the latter depend on N, so is the former. When 
hi converge as N —> oo, this dependence can be removed by using the continuity of p(x) in hi, 
established below in Proposition |5.8| 


Analogues of Theorem 5.3 are known in the literature, cf. |MMS| . [AGZj . jFdlj . Our proof of 
Theorem 5.3 relies on a different characterization for p maximizing the functional Iy[-\. 

The restriction 0 < p{x) < 0 _1 leads to the subdivision of M into three types of regions (we 
borrow the terminology from jBKMM j): 

• Maximal (with respect to inclusion) closed connected intervals where p(x) = 0 are called 
voids. 

• Maximal open connected intervals where 0 < p(x) < # _1 are called bands. 

• Maximal closed connected intervals where p(x) = 0~ 1 are called saturated regions. 

In a related context of random tilings and periodically-weighted dimers the voids and saturated 
regions are usually called frozen, while bands are liquid regions. 

We further define the effective potential Fy(x) through (cf. ©) 


Fy(x) = 20 


In | x-y\p(y)dy-V(x). 


Lemma 5.5. There exist k constants fi, ■ ■ ■, fk, such that 

(88) Fy(x) — f < 0, for all x in voids in [hi, 6/]; 

(89) Fy(x) — f > 0, for all x in saturated regions in [hi, 5/]; 

(90) Fy(x) — f = 0, for all x in bands in [hi, bi\. 

Proof. This characterization is readily obtained from the variation of the functional Iy\-\, cf. jSTj . 

|DS] . □ 

Further fix a parameter p > 2 and let p.^ denote the convolution of the empirical measure 
pjy with uniform measure on the interval [0,IV -p ]. The following statement is a generalization 
of Proposition 2.16[ Let us emphasize that the measure p here and below depends on the filling 
fractions hi = ni(N)/N, and thus changes with N. 
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Proposition 5.6. Let fi be the maximizer of Lemma 5.1. Then there exists C G 
all 7 > 0 and all N we have 


such that for 


P N (v(ii N ,n) > 7 ) < exp (CiV In 2 (AO - 7 2 iV 2 ), 


where V was defined in (35), (36). 

Proof. We start by analyzing the asymptotic of the formula for Pjv- We observe that the definition 
of Iv['] by (87) makes sense for discrete measures. It is important here that we integrate only 
over x 7 ^ y in (87) as otherwise the integral would be infinite. Using for the double product in 
the following corollary of the Stirling’s formula 

T( 'r(h,)' ) = h \ l + h-^+oo, 


(91) 


we can write for every ,... ,£n) G W 1 


N 


(92) F n (£ u ...,£ n ) = 

where 


exp 


(20N(N - l)ln{N) + N 2 I v [mes[h,... ,£ n ]\ +0(NIn 

Z N 


1 N 

mes[ 4 , ■ ■ -An\ = 


i =1 


N 


Let us obtain a lower bound for the partition function Zn in (92). For that let X{, i = 1,. 
be quantiles of (j, defined through 

i - 1/2 


,N 


fi(x)dx = 


N 


1 < i < N. 


Since /j,(x) <6 1 , 6{xi + \ — xf) > 1/N. Therefore, using the asymptotics of ai(N) and bi(N) of 
Assumption [T] we conclude that there exists an absolute constant U (independent of N) and a 
configuration (£ 1 ,... ,£n) £ W(y such that 

(93) \Nxi-ii\<U 

for all 1 < i < N except for 0(ln(N)) ones. Arguing as in Proposition |2. 16] we write 

jy 207 V(iv-i) - ex P (n 2 Iv [mes[ii,... ,l N ]] + O(NlnN)^ = exp (n 2 I v [h\ + 0(Nln 2 (N))y 

Note that the remainder in the last formula is 0(N In 2 (N)) instead of 0(N ln(N)) in the proof 
of Proposition 2.16 This is due to an additional error produced by the ti for which (93) does not 
hold: there are 0(ln(AT)) of such and each of them produces (at most) 0(N ln(A r )) of error. 


From here on the proof literally repeats that of Proposition 2.16 


□ 

Corollary 5.7. Recall the notations ||(/||i /2 and ||g||Li P of Corollary 2.12 Fix any p > 2 and let 
fi(x)dx be the maximizer of Lemma 5.1. Then there exists C?M such that for all 7 > 0, all N 
and all f we have 


(94) ¥ n 


g(x)n N (dx) - / g(x)ix(x)dx 


> 7115-111/2 + 


115-11 Lip \ 


NP J 


< exp 


( C Nln\N)-^pj. 




Proof. The proof is the same as for Corollary 2.17 


□ 


Proof of Theorem 5.3. The desired statement follows from Corollary 5.7 with 7 = q ■ N £ 1//2 , 
q > 0 . □ 
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Recall that the set of measures 0 and the equilibrium measure y of Lemma 5.1 implicitly 
depend on many parameters. For the next statement we reconstruct a part of this dependence in 
the notations. For h = (hi, ..., hk) let y n be the maximizer of Iy over 0 n , the set of probability 
measures y absolutely continuous with respect to Lebesgue measure, with density bounded above 
by 1/9 and with mass y([hi, 6,.]) = hi, Yli=i F, = 1- 

Proposition 5.8. Assume that 0 < hi < 9~ 1 (bi — hi) for all i = 1, ..., k and that the measure 
[i Ui has at least one band in each of the intervals [hi,bi], i = 1,... ,k. Then there exists a finite 
constant Ca such that for h! close enough to h 


(95) 


V(y n , y n ) < Cfi\\h — h'\ 


Remark 5.9. Arguing as in the proof of Corollary 2.17, one deduces from (95) that for any smooth 


function g(x), the averages f g(x)y n (x)dx and f g(x)y n ' (x)dx are close. Their difference goes to 


zero as n 


n. 


Remark 5.10. In all the applications that we present, the assumption that the measure y ni has 
at least one band in each of the intervals [hi, bi], i = 1 ,... ,k, is satisfied automatically, see the 
end of Section |9.1| for a general argument in this direction. 


Proof of Proposition |7Q[ By definition, for any probability measure u with density u(x) 

Iy(v) = I v (n n ) - 9V 2 (u,y n ) - J (v(x) - 29 J log \x - y\y n (y)dyj (v(x) - / x n (x))dx 

/ k 

Fy(x)(v(x) - n n (x))dx + ^2 fi(v - n n )([hi, S*]) 

»=l 

where Fy(x) = Fy(x) — fi on [hi, bi] is nonpositive on voids, nonnegative on saturated regions 


and vanishes on bands for y n , cf. (88)—(90). We next choose hi ^ h and let v be a probability 


measure with filling fractions h! so that 

Iv(v)<Iv(^')- 

The above decomposition and (v — y n ')([hi, bi]) = 0 implies 

ev\y‘,y) -jMvP'a) - Ax))* < A - /#K*>M*) - Ax))*. 

Observe that f Fy(x)(u(x) — y n (x))dx is nonpositive for any probability measure v whose density 
is at most 9~ x , and vanishes if u — fa n is supported on the bands of y n . We assume that this is 
the case to deduce that 

(96) V 2 (y W ,y h )<V 2 (u,y h ). 

Finally we choose u = yP + (L( — hi)^fdx where Bi denotes a subset in the interior of bands 

in [hi,bi\ so that on Bi we have 5 < dy n /dx < 1/9 — 5 for some small <5 > 0. Note that our 
assumption on the existence of bands in all the intervals [d/, bf implies that Bi is not empty. For 
n( — ni small enough (smaller than 6\Bi\), v e 0 n . Hence, we have 

07) a=- e [ f 

for a constant Ca > 0. Combining (96) and (|97|) we obtain 


In \t — s | dt ds < Cn\\h — n 


/112 
oo 


^(yf 1 , ii h ) < Cfi\\h — n , ||^ 0 . □ 
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5.2. Functions R and Q. We work with the equilibrium measure // of density p(x) in terms of 
its Stieltjes transform G(z) defined through 


(98) 


fC 

G»{z) = / 

J — C 


p(x)dx 


z — x 


Observe that (98) makes sense for all z outside the support of p(x) and G fl (z) is holomorphic 
there. 

We are going to extensively use the following two notations for any function F(z) of a complex 
variable z: 


F(z + iO) := lim F(z + i<J), 
6 — >o 


F(z — iO) := lim F(z — i<5). 
5->o 


At every point x where p(x) is continuous, it can be reconstructed by its Stieltjes transform via 


(99) 


h(x) = x-t (Gfj,(x + iO) - G^x - iO)), x G 
27T1 


On the other hand, the equilibrium measure characterization (90) implies (cf. [PS] , [ST]) that for 
x in a band of the equilibrium measure 


( 100 ) 


d 


0(G fl (x + iO) + G m (x - iO)) = — V(x). 


We also define the Stieltjes transform Gn(z) of the prelimit empirical measure ( 86 ) through 

/ °c 1 1 N ' J 

-^—ji N (dx) = — ^2 „ _ n /AT , (4, • ■ ■, 4v) e w % is FAT-distributed. 

-OO ^ ^ ^ 


The functions R^(z), Q^z) defined in (75), (76). are important for our asymptotic study. 
Proposition 5.11. Under Assumptions [71 111 [1 [ 5 | Rfi(z) is holomorphic in M.n- 


Proof of Proposition 5.11 . Since G^(z) is holomorphic everywhere outside the (real) support of 
the equilibrium measure /u(x), so is i? M (z). 

Further take any point x in the support of Choose a simple contour 7 in Ad at \ [&$,&*] 

enclosing x. Observe that as N —» 00 under the change of variables £ = Nz, we have 


( 102 ) 


1 - 


6 




= 11 - 1 - 


6 


1 + 


£ - 4 - 1 


N z-ii/N 

1 


= exp 


= 1 + T7- 


N z-li/N-l/N 


1 9 

N ' z - ii/N 

1 


= exp 


+o u 

9 


N z-ii/N 


+ 0{ N 


Now fix ni, i = 1,..., k and choose the filling fractions rn(N) in such a way that \ni(N) — Nhi\ < 1 
for all N and i. We claim that Rn(Nz) of Theorem |4 .1 1 converges to i? /t (z) uniformly on 7 . Indeed, 
the functions <frf] converge to (fr^ by Assumption^ By (102) and Theorem 5.3 the two expectations 
in the definition of R^(Nz ) ap proximate exp(— 9G u (z)) and exp(— 9G u (z)), respectively, where v 

with filling fractions ni ^ , ... . Since ni ^ —> hi as N —> 00 , 


5.1 


is the maximizer of Lemma 

Proposition 5.8 implies that G u (z) converges to G^(z 

Since uniform convergence on 7 of holomorphic (everywhere inside the contour) functions 
Rn(Nz) implies the same convergence inside the contour and the holomorphicity of the limit, 
we are done. □ 


Remark 5.12. Alternatively, one can prove Proposition 5.11 by showing that Rfi(z) is continuous 
near real x via combining (100) with (73). 
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Despite their similar form, the analytic properties of the function Q^z) are very different. 
Observe that 

(103) {Q»( z )) 2 = { R v( z )) 2 - ^ + (z)ct)~{z). 

Thus, Q^(z) is a branch of a two-valued analytic function which is the square-root of a function 
holomorphic in M.n- Note that Qn(z) is analytic outside Uf = 1 [oj, 6 j]. On the other hand, com¬ 
bining (100) with ( p73] ) one observes that Qn(z) must have discontinuities in the bands of /i(x) 
and the endpoints of the bands should be its branching points. Therefore, Q^{z) has at least two 
branching points inside each interval [a/, b,]. which also have to be zeros of Q^(z). For our method 
it is important to assume that there are precisely two zeros in each interval, as is summarized in 
Assumption [lj 


6. Second order expansion. 


The goal of this section is to prove a generalization of Theorem 2.5 in the setup of Section [3} 
Theorem 5.3 implies that Gn(z) — G /1 (z) vanishes as N —> oo uniformly over z in compact 
subsets of C \ (Ji=i[oj> bj\. Moreover, since all the involved functions are analytic in z, we can 
infer also the uniform convergence of the derivatives. 

Similarly to Section 2 . 2 , we introduce a deformed version of Gn(z). Take 2 m parameters t = 
(ti, •.., t m ), v = (ui,..., v m ) such that v a + t a - y / 0 for all a = 1 ,m and all y € uf = 1 [a*, 6 *], 
and let the deformed distribution P^ v be defined through 

N m 


(104) 


¥% v (£ u ...,£ n ) = 


Z N 


Z N (t,v) 


' Pjv(^1, ' • • ■ nn *+ 


i= 1 a=l 


Va-li/Nj ’ 


where Pat was defined in (70). In general, P^ v may be a complex-valued measure, but the normal¬ 
izing constant Z^(t. v) in (104) is chosen so that the total mass of P^ v is 1, i.e. mv \i) = 1 . 

The numbers t a , a = 1,..., k are always assumed to be in a small neighborhood of 0, which guar¬ 
antees, in particular, that that the measure is normaliza ble, i.e. Z/v(t, v) ^ 0 . 

Observe that P)y V satisfies the assumptions of Theorem 

§ 


4.1 


with (in the notations of Assumption 


«W = (^(^) + 4^(|)+o(T 


U N 


X = 


N' 


JY ) ( AT ) ^ ( AT2 


N 


N 2 


nb,-MA + ^ + D. 


a=l 


We also define 


Vw(a) = <Pn( z ) + Y1 


$ (z) 


Clearly, 


4>n( N z ) = + 


t a + v a - z 

^n ( z ) 


^n(z) = <Pn(z)+Y^ 


< 2=1 


4> + {z) 


V a - Z 


N 


+ 01 ^ 


E[ ( U “ “ Z ) ( ta + V a ~ z). 


We further define as the empirical distribution of P)^ v and set, cf. (18) 

N 


( 2=1 

3>t,V 


/ °° 1 1 w— 

^ z _ x »n ( dx ) = N J1 


1 


NZ^z-ii/N' 

1=1 ' 


{l \,..., In) is Pj^ v -distributed. 
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Define 

(106) 


A G N {z) = N(G N {z)-G tx {z)). 


We will now formulate a generalization of Theorem 2.5 For that we need to introduce certain 
notations from the theory of hyperelliptic integrals. 

Fix k simple positively-oriented complex contours 71 , 72 ,... , 7 ^, such that each 7 * (1 < i < k) 
encloses the segment [di,bi] (and thus, also [ct*,/?*]), but not the other segments. 

Let P(z) = po + piz + • • • + pk- 2 Z k ~ 2 be a polynomial of degree k — 2 , and consider the map 


n : P(z) HA ( ——7 


P(z)dz 


1 


P{z)dz 


2 ™ hi nil "' ’ 2?ri h nil v(*-<*)(*-w 

Note that the sum of the integrals in the definition of D equals (—27ri) times the residue of 

—7 - , P ^ = at infinity, which is zero. Therefore, H is a linear map between (k — 1)— 

rii=i v( 2 -«i)( 2 -^i) 

dimensional vector spaces. This map is very non-trivial, but it is known to be an isomorphism 
for any k > 2 (cf. [Du . Section 2 . 1 ]). 

Using 12 we can now define a more complicated map T. Given a (continuous) function f(z) 
defined on the contours 7 \ and such that the sum of its integrals over these contours is zero, we 
define a function T z [/] through 

ui/i = m + „ p(z) =, 
llt=l V( z -«*)(«- Pi) 

where P(z) is a unique polynomial of degree at most k — 2 , such that for each i = 1 , 

oTT / T*[/]dz = 0. 




27ri 


7i 


The polynomial P(z) can be evaluated in terms of the map fl via 


p = n 


-1 


( 2 jri £ -£i £ nz)dz ’ ""“£;£ /(2)<i £ ' 


'71 J 72 ■'7fc 

Note that the map / ha T z [/] is linear and does not depend on t, v. 

Theorem 6 . 1 . Fix m = 0,1,... and choose complex numbers v = (v\,... ,v m ) C (C \ 
Uj=i [df, bi]) m . Under Assumptions M-Lm we have as N —> 00 


(107) 


d r 


dt\ ■ ■ ■ dt m p v 


T, 


E p t,v (A Gn(u)) 

e- 1 


d m 

= o(l) + — - — 

t a =0,1 <a<m ut 1 • • • Ot r 


dz 


■ ^(z)e~ eG h) + ip+ (z)e 9G h) 


27rinf=i \/Ju - a i){ u - Pi) hLi'n ( u ~ z ) H h 

+ + <P + (z)e ea -^ (y - ») h“ lz) 


t a = 0, l<a<m 


where 7 $ are simple positively-oriented contours enclosing the segment [dj,£>j] (i/ie points u and 
vi ,..., v m are outside the contours). The remainder o(l) is uniform over u,v 1 ,... ,v m in compact 
subsets of the unbounded component of C \ U^ =1 7 j. 


Remark 6.2. Theorem 6.1 generalizes Theorem 2.5 in several directions. First, before we had 
k = 1 and now we allow any k = 1,2,.... Second, 9 was equal to 1 and now 9 > 0 is arbitrary. 
Finally, in Theorem 2.5 the functions H(z), cf ± (z ), G^(z) had a specific form. 
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Proof. Similarly to Theorem 2.5 the proof has two parts: algebraic manipulations and self¬ 


improving estimates for the remainders. The latter literally repeats Section 2.5 and we are not 
going to present it. 

For the former we start from the statement of Theorem |4.1| Making the change of variables 
£ = Nz, we can write for 2 outside U/=iRhj &i] 


(108) 


N 


n 1 


2=1 


£-4 


= exp 


N / 1 

E>»M 


^ 2—1 


9 


N z-ti/N 


= exp (-9G n (z) + ^ • £ G n (z ) + O (J^ 

= exp (-dG^z) - —A G n {z) + — • — G n {z) + 0 


Similarly, also 


N 


(109) n 1 + 


2=1 


£ - U - 1 


ex pfe ln ( 1 + AT 


S. 2=1 


= exp ^6Gfj,(z) + ^AG n {z) + 


N z-ti/N-l/N 

9 2 /2 -9 d 


N 


-G n {z) + 0\ w 


Recalling the definition of R il {z) 1 we conclude that the function R^(Nz) from Theorem 4.1 for 


jt,v 

■ N 


can be written in the following form 


(110) R n (Nz) = Y[(v a -z){t a + v a -z) ■ R[j{z) 


a =1 


+ (f) (z)e 9G ^ z) E p t,v fexp 


AG N (z)^J - 1^ +</> + (z)e 0G ^E F t, v ^exp (~^AG N (z) 


- 1 


, ^ N ( z ) -eajz) , ^n( z ) 6G n (z) 

+ N + N 

, (t>-{z)e~ eG ^ 9 2 d (j)+(z)e eG ^ (9 2 Q \ d 

+ N 2 d~z°^ Z) + N U~ 9 Wz°^ Z) 


+ °' N 


We further want to simplify the expression in the second line of (110), by replacing e h — 1 by h. 


under expectations. As in Theorem 2.5, for that we use the inequality 

\e h — h — 1| < |/i| 2 e^, /i£ C, 

and the fact that 


E F 


G n (z) 


exp 


G n (z) 


= ° ( jj ) » N^oo, 


which is established via self-improving estimates. 
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We therefore can rewrite (110) as 
(m) 

Rn(Nz) = Y\( v a-z)(t a +v a -z)- Rfj,(z)- ^^^ E ¥ t,y (&G N (z))+ ^ N }~\ ~ 6G ^ z) + ^ N } T ^ e 90 ^ 


a =1 


N 


cj) (z)e eG v( z ) O' 2 d cf) + {z)e eGli ^ f 9 2 

H-v;- ——(*fj,(z) + - 


N 


N 


N 


2 dz 


N 




+ ° [ N 


where the remainder is uniform over 2 in compact subsets of C \U-=1 [&iM- 

Let us now fix u outside the contours 71 ,..., 7 *,, divide ( TTTh by 27rin)ii(' y a — z)(t a + v a — z) ■ 
H(z ) • (u — z) (the definition of H(z) is given in Assumption |4l) and integrate over the union of 
contours 77 We get 


( 112 ) 


6 

2iri 


nil y/(z - ati)(z - Pi) 


Ui=i 7* 


u — z 


E p t,v (A Gn(z)) dz 

r N 


dz 


• r / ^TT,^■ ( ipZ(z)e~ eG ^ z) + 'ipt(z)e 9G ^ z) 

+ f(z)e-^W^G /1 (.) + ^(z)e^ {2} (y - d) ^G„(z) 

’-NR n (Nz 


+ 


NR u (z) . 

+ 7-^Avt-w + o(1) )dz. 


2iri Ju* =l7i \(u- z)H(z) (u - z)H(z) 


Since H(z) is non-zero on all the segments [a*, 6 j], by a proper choice of small contours 7 \ we can 
guarantee that H(z ) is non-zero inside the integration contours. Then the terms in the last line 
of ( 112 ) vanish as they have no singularities inside the integration contours. 

Turning to the first line of (112), note that E p t,v (A Gn(z)) is analytic outside the contours of 

integration and decays as 1 / z 2 when z -A 00 . Therefore, we can compute the integral as (minus) 
the sum of the residues at z = u and at 2 = 00 . The former is 


9 y/{u- a*) (it - Pi) ■ E p t,v (AG N (u )), 


2—1 


while the latter is a polynomial Pn(u ) of degree at most k — 2 (to see that one uses (z — it) 1 = 
z~ l Ylm>o( uz ) m )• We conclude that 


(113) 9 ■ Ept,v (A Gjy(u)) + 


Pn(u ) 


nil \/{u - CLi){u - Pi) 


dz 


, _ _ r , . TT . • \^{z)e~ eG ^ z) +'if}±(z)e 6G ^ z) 

2vri ■ nil \/{u ~ oci)(u - Pi) ~ Z ) H ( Z ) \ 

+ r(z)e- SG "^j^G„(z) + 4,+ (z)e eG -^ (A - s) A G(1 (A + 0 (l). 


Now we are in a position to apply the map T u . Indeed, the integral of Gjy(z) around 7 \ is 
deterministic and equals rii(N)/N. On the other hand, the integral of G^(z) around 7 * equals the 
total mass of /x(x) inside 7 j, which is h*. Since, h* = rii(N)/N, the integral of E p t,v (AGn(u)) 
around each loop 7 i vanishes. 
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Therefore, 

(114) E p t,v(A G n (u)) 
= T„ 


e 


-i 


dz 


2 vri • nil V( u - a i)( u - ft) AjLi 74 

# 2 d 


^(z)e~ eG ^ z) + if+(z)e 6G ^ 


+ fWe-^ (z) +^V G ^ } (V/2-0)|^(ft(ft) 


+ o(D. 


It remains to show that (114) can be differentiated with respect to ti ,..., t m and t\ = ■ ■ ■ = t m = 0, 


which is done in the same way as in Theorem 2.5 see Section 2.5 for the details. It is important 
to note here that due to its explicit definition, the map T u is continuous and even Lipshitz (in 
uniform norm on the contours 7 j, i = 1 ,,k). Moreover, it is linear and does not depend on t. 
Therefore, its appearance in the formulas does not affect the argument. □ 


7. Central Limit Theorem 


The following two theorems are corollaries of Theorem 6.1, cf. Section 2.2 


Theorem 7.1. Under Assumptionsthe joint moments of the random variables N(G]y(z) — 
^p n Gn(z)) approximate (uniformly in z in compact subsets of C \ Ui=i.[®j, h]) those of centered 
Gaussian random variables with covariance 

N 2 (e Pn (G n (z)G n (w)) -E Pn G n (z)E Pn G n (w)^ = ^ £(z,w) + o(l), N 00 , 


where 


(115) £(z,w) = T h 


Vll(*-<*)(*-ft) 


(w - ft 2 yn i( w - a i)( w - ft) V ( z ~ 


E 


+ 


2(z — w) \z — cti z — fi. 

Remark 7.2. The map T was defined in Section [ 6 } In cases k = 1 and k = 2 the resulting 
covariance function £(z,w) can be brought to a more explicit form. When k = 1, 


£( 2 , w) = — 


1 - 


zw - l(ai + Pi){z + w) + «i/3i 


(z-w) 2 l y '{z- ai)(z - (3i)yJ(w - ai)(w - ft) 


The expression for k = 2 involves the values of complete elliptic integrals and we do not pursue 
it here, cf. [ BDEl Section 3]. 

Remark 7.3. Note that the numbers ct/, /3/ may depend on p(x), which, in turn, may depend on N 
through filling fractions hi = ni(N)/N. However, this dependence is continuous (see Proposition 
5.8) and thus, if hi converges as N —>• 00 , then so is €(z,w ), and Theorem 7.1 turns into a 


conventional central limit theorem for NGn(z) as N —> 00 . In particular, if k = 1, then this is 
always the case. 


Remark 7.4. In the general beta random matrix models the covariance has precisely the same 
form, cf. |J3j, jSh^h) . |Bo(lu2| . 
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Corollary 7.5. Take m> 1 real functions on Uf = 1 [aj, 6 j] that can be extended 

to holomorphic functions in a complex neighborhood B o/uf = 1 [aj, bi]. Under Assumptions [f]-[5], as 
N —> oo the joint moments of the m random variables 

N 

C fj = J2(m) - v/j(^))> (^i) • ■ -An) is P N -distributed, 

Z— 1 

approximate those of centered Gaussian random variables with asymptotic covariance 
(116) Ep N C fi C fj = —/ / fi(u)fj{v)£(u,v)dudv + o( 1), 

\ Z7Tl ) «/Ui7i JUir/i 


where n, i = 1,..., k, are positively oriented contours in B that enclose [ai,bi\, respectively, and 
£(u,v) is given by (115). 


in the asymptotic expansion of Ep N f{if) for functions / holomorphic in a neighborhood of 
U i = i[ai,bi\. 


Remark 7.6. Similarly to Corollary |7.5|, Theorem 6.1 can be used to obtain the first two terms 


8. NON-VANISHING WEIGHTS 


Assumption [5] of the general setup of Section [3] was vanishing of the weight at the end points 
of the supporting intervals. For future applications it is convenient to relax this condition and 
replace it by the following exponential bound on probabilities of having particles at a* (IV) + 1 or 
bi(N)- 1. 

Assumption 6 . We require the existence of constants C\, C 2 , C 3 > 0 such that for all N = 
1 , 2 ,..., the PAr-probability of the event 

lj = ai(N ) + 1 or £j = h(N) — 1 for at least one pair 1 < j < N, 1 < i < k 

is bounded from above by Cd exp(—C^IV*- 73 ). 


Theorem 8.1. If we replace Assumption [5] by Assumption [ 6 | then the results of Theorem 6.1 
Theorem 7.1, Corollary 7.5 are still valid for measures Ptv- 


The proof of Theorem 8.1 is based on the following observation. 


Proposition 8.2. In the notations of Theorem f.l, suppose that the condition +1) = 

4>~j^(bi(N)) = 0, i = 1,. .., k, does not hold. If 4> N {f ) are holomorphic in a domain Mn C C, then 
so is R n (0 except for at most 2 k simple poles. These poles are at points {ai(N) + 1, bi{N) — 1 }^ =1 . 
Under Assumption [ 6 ] the corresponding residues decay exponentially in N as N —>• 00 in the same 
sense as in the bound of Assumption^ (perhaps, with different constants C\, C 2 , C 3 ). 

Proof. We repeat the proof of Theorem |4.1| and observe the same cancelation of the poles. The 


only poles for which the cancelations do not occur, are endpoints of the interval: in Theorem 4.1 
the functions (j)± were vanishing at these endpoints, but this is no longer the case. The residue 


at an end-point can be bounded for the first term in (78) by the probability to have a particle at 
such an end-point (denote it by m ) multiplied by 


max 


n 

z: £i^m 


1 - 


6 


m — li 


< 


N -1 


1=1 


n ( i+ « 


= N. 


Thus, the exponential decay of the probability in Assumption [ 6 ] implies the exponential decay of 
the residue. For the second term in (78) the argument is the same. □ 
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Proof of Theorem 8.1. Note that an analogue of Proposition |8.2| is readily established also for the 


deformed measures as in SectionpJ Indeed, we need such measures only for small (i.e. tending 
to 0) values of t/, but then the exponential bounds on probability remain valid. Therefore, all 
the arguments of Section[ 6 ]go through for the measures Pat. Indeed, the only difference between 
Theorem 4.1 and Proposition |8.2| is in the appearance of finitely many simple poles. However, 
since the residues are exponentially small, these poles will only add exponentially small terms 
to all the remainders and thus will not contribute to the expansions in powers of 1 /N that we 
study. □ 


9. Examples 


The aim of this section is to demonstrate how the Assumptions [l}|4] are checked in applications, 

Theorem \7.1\ and Corollary |7.5| for certain stochastic 


which yields the validity of Theorem |6.1 
systems. 


9.1. Multi-cut general 9 extension of the Krawtchouk ensemble. The first example is an 
extension of that of Section [2] to general values of 9 and k. 

We fix k = 1, 2, ... and take 3 k numbers a/, 6 /, hi , such that 

a\ < bi < a.2 < l>2 < ■ ■ ■ < hk < bk, 

0 < hi < 9~ 1 (bi — hi) for all i = and Yli=i = T Then we further choose numbers 

a,i(N), bi(N), rii(N) such that the model fits into the setup of Section [3} The weight w(x]N) is 
then defined through the identity 


(117) 


w(x; N) 
w(x — 1; N) 


t-t x — b{(N) 

-j-1 x + 1 - di(N )' 
1=1 


Note that (117) agrees with the conditions w(x-,N ) > 0, x € ujL 1 [a/(A r ) + 1 ,bi(N) — 1] and 
with Assumption pi Since the weight w(x;N) is supported on U* = 1 (aj(lV), &j(JV)), we need to 


supplement (117) by the choice of k constants 
(118) 


(AT) = w(ai(N) + 1; N), i = l,...,k. 


Observe that the multiplication of all q by a same constant leaves the probability distribution 
unchanged. Therefore, for each N the system depends on the choice of 4/c — 2 constants. In 
particular, if k = 1 and 9 = 1, then (up to a shift of the lattice) we arrive at the example of 
Section [2j 

If we now assume that all the parameters are chosen so that as —> oo 

ai{N) = Nhi + 0{ 1), bi(N) = Nbi + 0( 1), m(N) = Nhi + 0(1), 


Cj(JV) = exp(ATci), i = l,...,k, 

then using Stirling’s formula for the factorials appearing in the explicit expressions for w(x;N), 
it is easy to see that the model satisfies Assumptions [lj [2j [3j Further, 

k k 

(/) + (z) = ~Y[( Z -bi), fr{z) = X\(z-ai)- 

1=1 1=1 

The function Rfj,(z) is an analytic function in z, which is 0(z k ^ 1 ) as z —> oo. Therefore, by 
Liouville’s theorem Rfi(z) is a polynomial of degree at most k — 1. Hence, the quadratic equation 
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(103) implies that Q^{z) is the square root of a degree 2k polynomial. In other words, decomposing 
into linear factors we get 


(119) 


k 


Q„{z) = 2 


N 


IK*- 


&i)(z 


Pi). 


On the other hand, our choice of hi guarantees, that the equilibrium measure /r( x)dx in each 
interval [hi,bi] has some points where fi(x) > 0 and some points where fi(x) < 0 _1 . Combining 
this with the observation (that is immediate from (|76[) ) that Q^(a*) and Q^(bi) have different 
signs for each i = 1 we conclude that a*,/?* E [hi,bi\. For each i there are two options: 

either fj,(x) has a band inside [cn, bi ], and then an and /3j are endpoints inside this band, or [hi, bi\ 
is a union of a void and a saturated region. Let us explain that the latter is impossible. Indeed, 
then there must be a point x E M and e > 0 such that /i(x) = 0 on (x — e, x) and n(x) = 6 on 
(x,x + e) (or vice versa which is considered in the same way). Therefore, as z approaches x along 
the real axis from the left, the Stieltjes transform G^[z) explodes: 

lim GHz) = lim [ ^ ^ dt = — oo. 

z — r Z—>■£— I Z — t 


This implies lim z _>. 3 ;_ exp(G M (;z)) = 0 and lim z _>. x _ exp(— Gp(z)) = +oo. But then the dehnition 
f75| ) of Rfj,(z) implies that R tJ {z) has a singularity at z = x, which contradicts the fact that 
Rn(z) is a polynomial. Note that it is crucial in this argument that ((H does not have zeros inside 


(J(a,;, bi), and indeed the end-points hi, bi might separate voids and saturated regions. 

The conclusion is that this class of probability models satisfies Assumptions [l]|5] and Theorem 
6.1, Theorem 7.1, Corollary 7.5 are valid for them. As far as we know, these results are new with 


the exception of the case 6 = k = 1 . 


9.2. Lozenge tilings. The second example demonstrates that at least two instances of proba¬ 
bility measures arising in the study of uniformly random lozenge tilings fit into our framework. 

Consider an A x B x C hexagon drawn on the regular triangular lattice. We tile this hexagon 
with three types of elementary lozenges (which are unions of adjacent triangular faces of the 
lattice), cf. Figure [3j There are finitely many such tilings, and we are interested in the asymptotic 
behavior of uniformly random tiling as A,B,C —> oo. This is a well-studied model, with many 
results available, cf. [CLP] . [BKMM] . [JN] . [G], [BoGoj . [FT], [P2] . 

Let us dissect the hexagon by a vertical line at distance t from the left boundary. There will be 
a fixed (depending on A, B, C, t ) number N of horizontal lozenges on this line; let Pat denote the 
distribution of these lozenges. This distribution can be computed by noticing that the tiling can 
be viewed as two Gelfand-Tsetlin patterns glued together, as shown in the right panel of Figure 
[3] The enumeration of Gelfand-Tsetlin patterns is well-known and can be used to compute the 
distribution Fjy (see [CLP ]. [G], [BP] . [BuGo] for more details). Assuming t > max(£>, C) as 
in Figure [3] which yields N = B + C — t, and introducing the coordinate system such that the 
lowest possible position for horizontal lozenges on the tth vertical line is 1 and the highest one is 
A + B + C — t, we obtain the formula 

(■A + B-|-C , + 1 — t — ii)t—B i)t—C j 


(120) Piv^,... M=^ n& - £ ;) 2 n 


where (o) n is the Pochhammer symbol, (a) n = a(a + 1) • • • (a + n — 1), and Zjy is a normalizing 
constant (which can be computed explicitly in this case). The distribution of the form (120) is 
known as Hahn orthogonal polynomial ensemble. 
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Figure 3. Left panel: Lozenge tiling of the 3x4x5 hexagon and 3 horizontal 
lozenges on the sixth from the left vertical line. Right panel: two Gelfand-Tsetlin 
patterns (inside blue and red contours) corresponding to each tiling. 


Take a large parameter L and suppose that 
(121) A = AL + 0( 1), B = BL + 0{ 1), C = CL+ 0(1), t = tL + 0{ 1), 
for positive constants A, B , C, t. We assume that t > max(L>, C) — other possibilities for t 


are considered similarly. Let us check that under such choice of parameters the ensemble (120) 
satisfies Assumptions [T] - [5} 

• Assumption [l] is satisfied due to Stirling’s formula applied to Pochhammer symbols. The 
potential V(u) has the form: 

V(u) = — {A + C — u) ln(A + C — u) + (A + B + C — t — u) ln(A + B + C — t — u) 

— (t — C + u) ln(t — C + u) + u ln('u). 

• Assumption ^1 is empty, as there is only one filling fraction, ni(lV) = N. 

• Assumption [3] is immediate from the definitions, and we have 

4 > + (z) = (t — C + z)(A + B + C — t — z), <t>~(z) = z(A + C — z). 

For Assumption [ 4 ] note that B 1i (z) is an analytic function which grows as 0{z 2 ) as z —> 00 


and therefore it is a polynomial of degree at most two. Hence, by ( |103 ), the function Q^{z) 
is a square root of a polynomial. The definition implies that Q^{z) is O(z) as z -> 00 . 
Thus Qn(z) is a square root of degree two polynomial and has the form 

Qn(z ) = const • \J{z — a)(z — b). 

The points a and b are necessarily endpoints of the band of the equilibrium measure 
n(x)dx. 

Assumption [ 5 ] is immediate from the formula (120). 


The conclusion is that Theorem 6.1 Theorem 7.1 Corollary 7.5 are valid. This implies a form 
of the Central Limit Theorem for fluctuations of lozenge tilings of a hexagon. We remark that 
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Figure 4. Lozenge tiling of the 4x7x5 hexagon with a rhombic 2x2 hole (shown 
in blue). The center of the hole is at distance t from the left side of the hexagon. 
The remaining horizontal lozenges on tth vertical line are shown in gray. 


the same CLT (and even stronger statement concerning joint asymptotic Gaussianity for several 
values of t) can be also established by other methods, cf. |P2] . |BDj . |BnGo2j . 


Remark 9.1. One can probably use the fact that R^(z) is a degree two polynomial to find an 
explicit formula for R^z) and thus, also for G^(z), and for the equilibrium measure / i(x)dx 
describing the limit shape for lozenge tilings, cf. Section 9.4 for a similar argument. Explicit 
formulas for p(x) were previously found by other methods in [CLP , [BKMM], |E3, m , EH- 


One can also analyze tilings of more complicated domains. Let us cut a rhombic D x D hole in 
the hexagon, as shown in Figure [4j Assume that the bottom point of the hole is at distance t from 
the left side of the hexagon and at height H (counted from the bottom of the hexagon along tth 
vertical line). Let Pw be the probability distribution of the horizontal lozenges (outside the hole) 
on tth vertical line induced by the uniform measure on all tilings of the hexagon with the hole. 
We can repeat the same argument as the one used for the complete hexagon for computing P^v- 
Assuming t > max(L>, C ) as in Figure [3j which yields N = B + C — D — t, and introducing the 
coordinate system such that the lowest possible position for horizontal lozenge on the tth vertical 
line is 1 and the highest one is A + B + C — t, we obtain 
( 122 ) 


N 


Ptv(G, • • • An) = 


Z N 


11 




i— 1 L 


(A + B + C+l-t-ejt-B Ai)t-c (H-li) D ( H-li) D 


We also need two filling fractions m and n- 2 - we consider only such tilings that there are m 
horizontal lozenges (on tth vertical line) below the hole and 77-2 lozenges above. Take a large 
parameter L and suppose that in addition to (121) we have 


H = HL + 0(1), D = DL + 0(1), H > 0, D> 0, H + D < A + B + C-t, 
ni = hi-L + 0(1), ri2 = fi2L + 0(1), 0 < h\ < H, 0 < ?i2 <A + B + C — t — H — D. 
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We again assume that t > max(B, C ) and remark that other possibilities for t are considered sim¬ 
ilarly. Let us check that under such choice of parameters the ensemble ( 122 ) satisfies Assumptions 

El-El 

• Assumption 1 is satisfied due to Stirling’s formula applied to Pochhammer symbols. 

• Assumption 2 is satisfied due to restrictions on hi, h 2 . 

• Assumption 3 is immediate from the definitions, and we have 

(j)+( z ) = (f - c + z)(A + B + C - t - z)(H - z ) 2 , </T{z) = z(A + C - z)(H + D - z) 2 . 

• For Assumption [i] note that R^z) is an analytic function which grows as 0(z A ) as z -A 00 
and therefore it is a polynomial of degree at most four. Hence, by (103), the function 
Qn(z) is a square root of a polynomial. The definition implies that Q^(z ) is 0(z 3 ) as 
z —> 00 . Thus, Qn(z) is a square root of a degree six polynomial. Further observe that 
according to definitions Q^{ 0) < 0, Q^(H) > 0, Qn(H + D) < 0, Qn(A + B + C — t) > 0. 
Together with analyticity of Qn(z) outside the support of the equilibrium measure this 
yields 

Qn(z) = const (z - c)\/(z- ai)(z - bi)(z - a 2 )(z - b - 2 ), 

where 

0<ai<bi<H<c<H + D<ci 2 <b 2 <A + B + C — t. 

As before, the points a\, 61 , 02,62 can be identified with the endpoints of the bands of the 
equilibrium measure. 

• Assumption [ 5 ] is immediate from the formula (122). 


The conclusion is that Theorem 6.1 Theorem \7.1\ Corollary |7.5| are valid. This implies a form of 
the Central Limit Theorem for fluctuations of lozenge tilings of a hexagon with a hole (with fixed 
filling fractions above and below the hole). As far as we know, these results are new; the same 
applies to the examples of the next sections. 

9.3. Arbitrary convex potential on M with no saturation. For our third example take a 
real convex analytic function V(x), i.e. such that V"(x) > 0 for all i£l. Fix a constant k > 0 
such that 

kV{x) 


(123) 

and consider a probability distribution 


lim inf 
x^roc 20m\x\ 


> 1 


(124) Pjv(4,- 

on A-tuples i\ < < 


Jn) = y 

ZjN 


n 


r (£j -£i + i)r (ij -u + 6) 

r (tj - - u +1 - e) J 


N 


ri exp 


2=1 




A i E Z. 


1<2<j<AT 
< £ tv such that 

£i = A* + Oi, Ai < A 2 < ■ ■ ■ < Ajvj 

In other words, we are in the framework of Section [3] except that now the distribution is sup¬ 
ported on an infinite and unbounded subset of M. N . In particular, when V(x) = x 2 , we obtain a 
discretization of the celebrated general (3 Gaussian ensemble from random matrix theory. 


Proposition 9.2. For any real analytic V{x) satisfying V"(x) > 0 and (123), the functional I K y 
of (87) has a unique maximizer (equilibrium measure) on the the set of all absolutely continuous 
probability measures on M (without any restrictions on the density). The equilibrium measure is 
compactly supported, has a continuous density fj,{x) and has a single band. At the end-points of 
the band the density of the equilibrium measure behaves like Cyfx. 

Proof. This is well-known, see m, HQ and references therein. □ 
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Note that the equilibrium measure /j,(x)dx of Proposition 9.2 is unchanged when we multiply 
both 0 and k on the same positive constant (as the functional is then multiplied by the same 
constant) and, therefore, depends only on their ratio k/6 . In particular, if we fix ratio and then 
choose small n (equivalently, small 9), then the density fi(x) will be smaller than 9~ l for all 0. 
But then the solution to constrained maximization problem with 0 < n(x) < 9~ l as in Theorem 
5.3 will be the same as the solution of the unconstrained minimization problem as in Proposition 
9.2 The conclusion is that for small n the equilibrium measure has no saturated regions; this is 


crucial for our considerations. 


We now want to show that Theorem 6.1 Theorem 7.1 and Corollary |7.5| apply in this situation. 


This is done by localizing the probability measure Pat onto a finite interval with help of the 
following proposition. 


Proposition 9.3. Assume that V(x) is convex and satisfies (123). Then with exponentially high 
probability the measure (124) is supported on configurations in a linearly growing interval, i.e. 


there exist constants C,D > 0 such that 


P n (-D< — < — 
1 ~ N ~ N 


< D > 1 — 


1 

C 


exp(— NC), 


N = 1 , 2 ,.... 


Proof. This a particular case of Theorem 10. 1| We remark that statements of this flavor are 
well-known as large deviations principles for the largest/smallest particles in both discrete and 
continuous log-gases, cf. pa Theorem 2.2], [Fe, Theorem 4.2], A( !Z . Section 2.6]. □ 


Proposition 9.3 motivates the definition of measure Pat as Pat conditioned on the event that 


\£i/N | < D + 1 for all i = 1,..., N. An advantage of the measure Pat is that ii/N are PAr-almost 
surely bounded, and therefore, there exists a finite complex contour which enclose them all and 
we can apply the developments of previous sections. 


Lemma 9.4. Assume (123) and that analytic V(x) satisfies V"(x) > 0 for all x E M. Then the 


results of Theorem 


6.1 


Theorem 


7.1 


Corollary 7.5 are valid for measures Pat. Here k = 1, the 


points a\ and b\ are two endpoints of the band of the equilibrium measure, 


0) = 1, T n {z) = 0, (f + (z) = exp(— kV'( z)), = 


kV"(z) 


exp(— kV'(z)). 


Proof. We will apply Theorem |8.1[ and we need to check all its assumptions. 
With w(x;N) = exp(—KNV(x/N)), we have 


(125) 


w(x ; N) 
w(x — 1; N ) 


- = exp 


kN V 


e-i 

N 


- V 


Therefore, the formulas for and tp + are obtained through large N expansion of the right-hand 
side in (125). 

We only check that Assumption [4] is satisfied, as the rest is automatic. We have 
Qn(z) = exp (- 6 Gfj,(z)) - exp(-KV , (z))exp(9G fl (z)). 


As before, Q^z) is a square root of an analytic function due to (103). Further, Q^ is analytic 
outside the single band of the equilibrium measure p,(x)dx. Near the end-points of the band, the 
density n{x) behaves like a square root and thus similar behavior for Q^i^z). Therefore, 

Qu( z ) = H ( Z )V( Z ~ a )( z ~ h )i 


where a < b are endpoints of the band. We will now show that H(z) has no zeros on M, which 
would also imply that H(z ) is holomorphic, since H 2 (z) is. 

We already ruled out H(a) = 0 and H(b) = 0, and there are 3 more cases to consider: 
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H(x) = 0, a < x < b. Then also Q^{x + iO) = Q^(x — iO) = 0. Therefore, 
exp(nV' (x)) = exp(2 9G^(x + iO)) = exp(20G ^{x — iO)). 


Comparing with (100) (V'(x) should be replaced by kV'(x) to match the notations) we 
conclude that 


exp 


(oG^fx + iO) - QG,jXx - i0)^ = 1. 


But this contradicts 0 < p(x) < 6 1 and (99). 


(126) 


For x > b, if H(x) = 0, then also Q(x) = 0 and 

exp(nV' (x)) = exp(2 OG^x)). 


But note that ( 126[) h olds at x = b as a corollary of (100). And as x > b grows, the 
left-hand side of (126) also grows due to V"(x) > 0, but the right-hand side decays due 
to the definition of G^{x). Thus, (126) can not hold at x > b. 

• For x < a the argument is similar. □ 

As a corollary we obtain a central limit theorem for the original measure Pat. 


Corollary 9.5. Assume (123), V"(x) > 0 for all and consider the probability measure P/v 

given by (124). Suppose that k is so small that the equilibrium measure has one band. Take any 
m functions fi(z ),..., fm{z) on M which are bounded and extend to holomorphic functions in a 
complex neighborhood of (—D,D). Then the m random variables 

N 

= ^2 {fj(ti) - Ep N fj{Zi)) , (h ,...,£n) is V N -distributed, 

i=1 

converge (in distribution and in the sense of moments) to centered Gaussian random variables 
with covariance 


(127) 


Jim 


P f fi(u)fj(v)C(u,v)dudv , 

7 ^7 


where C(u,v) is given by ( 22 ) with a± being endpoints of the band of the equilibrium measure 
p,(x)dx, and 7 is a positively oriented contour which encloses [a_,a+]. 


Proof. Due to Proposition 9.3 and boundedness of functions fj, the joint moments of J2iLi fj(^ 


with respect to P,y a nd with respect to P/v differ by exponentially small (in N) error. It remains 
to use Corollary 


7.5 


for P/v- 


□ 


Remark 9.6. One way to construct functions fj satisfying conditions of Corollary 9.5 is to take 
any analytic functions on \—D,D\ and extend them to all of 
[—D,D\. The condition that fj must be bounded on 


by setting equal to 0 outside 

can be weakened. However, we need 

N 


a growth condition on these functions, as we want to be able to replace Pp w (^/=i by 


fj(^i)Y w hb negligible error. The analyticity assumption can also likely be weakened, 
but we do not address this in the present paper. 


Remark 9.7. The covariance (127) has the same form as for random matrices and log-gases in 


the one cut regime. It depends only on the restrictions of functions fj onto the interval [a_,a+] 
and can be rewritten in several other equivalent forms, cf. HU Theorem 4.2], PS Chapter 3]. 
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9.4. (z, w)— measures. Our last example originates in the asymptotic representation theory of 
unitary groups U(N), cf. [O], jBOlj . [02| . 

Fix two sequences of non-real parameters z (IV) and w(JV) and dehne 


(128) P n (£i,...,£n) = ^- 

Zjn 


n 


T(tj -£i + l)r (£j -£i + 0) 

r (Ej - -li + l-6) 


N 

n 


^ r(z(JV) - £i)T(z(N) - £i)T(w(N) + £i)T{w(N) + if) 

on iV-tuples < £2 < • • • < In such that 

i% = Aj + Qi , Ai < A2 < ■ ■ ■ < Ajv ; Aj E Z. 

Here z(IV) and w(lV) are complex conjugates of z(IV) and w(IV), respectively. 

If 9 = 1, z(lV) = z + N + 1, w(IV) = w with Re(z + in) > —1/2, and {Aj }(^ 1 are identified with 
highest weights of irreducible representations of U(N), then (128) describes the decomposition 
of the character of the “generalized bi-regular” representation of the infinite-dimensional unitary 
group U ( 00 ), see 0 - 
With the notation 

w[x\ N ) = 


we have 

w(x; N) 


r(z(fV) - s)r(z(JV) - x)T(w(N) + x)T(w(N) + x ) ’ 
(x — z(N))(x — z (N)) 


w(x — 1; N) (x + w (N) — l)(x + w(IV) — 1) 


= 1 - 


w (N) + w (IV) + z(N) + z (IV) — 2 


+ oh 


x z 


The last formula implies that as \x\ —>■ 00 , the weight decays as 
(129) w(x; N) = o( | x |w(iV)+w(V)+z(A)+z(iv )-2 


Therefore, the real part of w (N) + w (N) + z(N ) + z (N) needs to be large in order to guarantee 
that the measure Pw is finite. Let us assume that as N —> 00 


(130) 


w(N) = Wqo ■ N + 0(1), z(N) = Zoo ■ N + 0(1), Re(woo + z^) > 1. 


Then w(x\ N) decays fast enough so that a condition of the form (123) is satisfied. Note that 
the original representation-theoretic case z(IV) = z + N + 1, w(fV) = w does not satisfy this 
assumption, and indeed the equilibrium measure in this case is known to be somewhat degenerate, 
and one generally does not expect to see the Gaussian behavior. This case with 6 = 1 was studied 
in [HOT] , 


Proposition 9.8. Assume (130). Then with exponentially high probability the measure (128) is 
supported on configurations in a linearly growing interval, i.e. there exist constants C,D > 0 such 
that 

r N (-D<%<^<D]> 1 - i .expHVC), 


N ~ N 


C 


N = 1,2,. 


Proof. This is a particular case of Theorem 10.1 


□ 

Proposition |9.8| implies that we can use for P^v the techniques developed in Sections [3]-[6j Let 
us present the functions and Q^, as they can be found explicitly in this case. Indeed, we have 

R u(0 = (£ - z oo)(£ - Zoo) exp(0G m (£)) + (£ + Woo)(£ + woo) exp(— 0 G M (£))- 
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Since R fl (£) is analytic and grows as 2£ 2 as |£| —> oo, is a degree two polynomial, i.e. 

Rh(Z) = 2£ 2 + AZ + B. 

Let us find the coefficients A and B. Expand in power series near £ = oo as 

= L + ^2 + 0(£ 3 ), 

where pi is unknown. Plugging this into the definition of i?^(£) and expanding up to 0(£ _1 ) we 
get 


(131) Rn(Q = (£ - Zoo )(£ - Zoo) 1 + 7 + 


6 0p\ + Q r - 


e 


+ 


(£ + Woo)(£ + Woo) ^1 - | ^ + 0(£ x ) 

= 2£ 2 + (Woo + Woo - Zoo - Zoo)£ + ZooZoo + W^Woo - 0(Zoo + Zqo + Woo + Woo) + 2 6 2 . 


Therefore, we can also find Q^ through (103): 

(132) Q^) = \J ( R M (£)) 2 - 4(£ - Zoo)(£ - Zoo)(£ + Woo)(£ + w^) = c^J - a_)(£ - a+) 


with explicit a± and c which are found by plugging (131) in (103) (the resulting formulas are 
somewhat complicated and we omit them). This is precisely the form we need for Assumption |4j 
Now repeating the argument of Lemma |9.4[, we see that the results of Theorem | 6 .H Theorem 


7.1, Corollary 7.5 are valid for measures P,v conditioned on the event of Proposition |9.8[ Thus, 


as in Corollary 9.5, we arrive at the following result. 


Corollary 9.9. Assume (130) and consider the probability measure P/v given by (128). Take 


m> 1 functions fi(z ),..., f m (z) on M that are bounded and extend to holomorphic functions in a 
complex neighborhood of (-D, D), where D is given by Proposition [5^. Then m random variables 

N 

C fj = ^2 - E P N fj(£i)) , (h ,...,£n) is P N -distributed, 


i =1 


converge (in distribution and in the sense of moments) to centered Gaussian random variables 
with covariance 

(133) lirn Ep N C fi C fj = 


N—>oo 


ri)^// f i ( u ')h( v> > C ( u ’ v ') dudv ’ 


(2tt1) J'y J 'y 

where C(u,v) is given by ( | 22 | ) with a± found from ( |131 ), ( 132| ) , and 7 is a positively oriented 
contour which encloses [a_,a+]. 


10. Exponential bound on the support. 


Take a continuous function V : M — > M and numbers e > 0, H > 1 + 9 1 such that 

(134) Yj' X ) 1 > 1 + e, when |x| > H. 

26 m |x| 

Set T = \6\ + 1 and assume that V(x) is increasing for x > H, decreasing for x < — H and is 
Lipshitz with a constant s for \x\ < H + T. 

Consider a probability distribution 


(135) P/v(4, • • •, £n) = 

Zjn 


n 


V(£j -£i + l)V(£j -£ i + 9)yr 

T{£j - £i)T(£j -£i + l- 6 ) 11 


exp 


~ N ' V ^ 
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on IV-tuples l\ < li < ■ ■ ■ < In such that 

(136) £i = Xi + 9i, Ai < A 2 < ••• < Ajv, A; G Z. 

The aim of this section is to prove the following statement describing the tails of Pa- 

Theorem 10.1. There exist two constants C = C(Q,s, H,s) and D = D(0,e, H,s) that depend 
only on 0,e,H,s, and such that 


Pat, n [-D< ( ^<^<d)> 1 - ^ • exp(-iVC), 


IV = 1,2,, 


The proof of Theorem 10.1 borrows ideas from similar proofs in |.I3j . |AGZI , Section 2.7], |Fej . 
but additional care is required because of the shifts by 6 in the definition of li. We present the 
proof as a series of lemmas. 

It is useful to consider several modifications of the measure Pat, which we now introduce. The 
probability measure Pa + is defined on the same space of Is by the formula 
(137) 


Piv,+ (fu, • ■ • An) = 


1 


n 


T{ij - u + i)r(<, - ti + e) 


N 


z «,+ isihjv m - wt, a 


i=l 


]~J exp ( —N ■ V 


N + l 


The probability measure Pa,++ is defined on the same space of & by the formula 
(138) Pjv,++(fu, • ■ ■ An) 




N 


Zn , 


++ 




i=1 


N + l 


lattice, namely 
(139) li = 


A i + Qi, 


p(k) 

p(^) 

rn(fc) 

r N > 

r N,+ 

> 

p(^) pW 

m(h) 

A ’ 

■ A,+ > 

V,++’ 


i < 

k, 

+ 1)) 

i > 

k, 


Ai < A 2 < • • • < A tv, A i G Z. 


We also let Z™, zf Z^ ++ to be the normalizing constant for the corresponding measures. 
Further denote 

Mat(Tl, • • • An) = ZjvPjv(4, ■ • • An), 

and similarly for Pa,+ , Pa,++, P^, P^ + , P^ ++ . 

All the constants ci,c 2 ,... in the following statements depend only on 6 ,e,H,s, the exact 
values of the constants might change from statement to statement. 

Lemma 10.2. Then there exists c\ > 0 such that for any 1 < k < N + 1, and any l satisfying 
( |136 ) there exists l' satisfying ( 139| ) , for which 

(140) |li - l[\ < T, for all l < i < N 

and 

M %\?) 


(141) 


< ci exp(Nci). 


M N (l) 

There also exists a (possibly different) l" such that 

(142) \li-l'l\<T, for alll<i<N 
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and 


(143) 


■S’ 


(£" 


M n (£) 


> — exp(— Nci). 
ci 


Similarly, for any £! satisfying (139), there exists i satisfying (136) such that (140) and (141) 
Further, for any i" satisfying (139), there exists i satisfying (136) such that (142) 

u(k) 


hold. 


and 


ik) 

(143) hold. Finally, the same statements hold for the measures P)y' + and P 


N,++- 


Proof. We will only prove the first two statements of the Lemma, as the rest can be proven 
similarly. 

Let x = (xi,..., xn) £ and y = (yi,..., yw) € be such that x 3 — Xj> 6{j — i), for 
1 IS i < 3 < N and there exists m = 1,..., IV and Mgl such that 


Vi = 


Xi, 

Xi + M, 


1 < i < m, 
m < i < N. 


We claim that there exists C 2 = C 2 (M) such that 


(144) —exp(—IVC 2 ) 

C2 


s n 

1 


r (Xj -Xj + l)r(xj - Xi + 9) F(yj - yi)T(yj -yj + 1-9) 
r (Xj - Xi)T(xj -Xi + 1-9) V(y 3 - y t + l)T(y 3 - y% + 0) 


< c 2 exp(lVc 2 ). 


Indeed, (J91J) implies that 

n 


r (Xj - Xi + l)T{xj - Xi + 0 ) T(yj - yi)T(yj - y* + 1 - 6) 


n 1+0 


Xj — Xj 


t<j - Xi)T{ Xj -Xi + l-6) T{ yj - yi + l)r(y, - Vi + 6) 7 ,< m< . 

Since Xj — Xj > 9 {j — i ), the last product is bounded below by ^exp(—IVC2) and above by 
C2exp(—IVC2) for some C2 > 0 . 

Take £ satisfying (136). Suppose that £j < —R for i = 1,..., mi, \U\ < H for i = mi+1,..., m 2 , 
and tj > H for i = m 2 + 1,..., N. 

Define £ = (£\,..., In) through 



£^1 

II 

i < k, 


[£i + 9, 

i > k. 

And further define £' = (£\,.. 

., £' n ) through 



1 

II 

i < m 1 


1 ^2? 

i > m 1 


We claim that (141) holds. Indeed, the ratio of the factors in double product Hi<j is bounded 


by two applications of (144), and it remains to bound 

N 


]Q exp ( NV 


i =1 


- NV 


N 


If i < mi, then £[ < £j and the monotonicity of V(x) implies that corresponding factors are less 
than 1. If * > m 2 , then £\ > £j and again the monotonicity of V{x) implies that corresponding 
factors are less than 1. Finally, if mi < i < m 2 , then the Lipshitz property of V(x) gives the 
desired bound. 
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Next, we construct £". For i < m\ we set £'[ := £ t . For i > m 2 we set £'■ = £i — T. Since 
H > 1 + we can always choose the remaining coordinates i'f, m\ < i < m 2 in such a way 
that l" satisfies (139) and (140). We claim that (143) holds. Indeed, the ratio of the factors in 
double product II i < ■ is bounded by (144), and it remains to bound 


N 


]Q exp ( NV 


i =1 


- NV 


If i < mi , then £! i > £i and the monotonicity of V(x) implies that corresponding factors are greater 
than 1. If i > m 2 , then £! i < ii and again the monotonicity of V{x) implies that corresponding 
factors are greater than 1. Finally, if rri\ < i < m 2 , then the Lipshitz property of V(x) gives the 
desired bound. □ 

Lemma 10.3. There exists C2 > 0 such that 


(145) 


Z 


(fc) 

N,+ 


Zn,-\ 


< C 2 exp(c 2 lV), l<k<N + l. 


Proof. For £' satisfying ( |139 ) let s(£) denote the corresponding £ of Lemma 
4,+ = E m a?+( 0 ^ c i wWa) E Mjv ;+ (s(/)) 


10.2 


Then 


< 


ciexp(ATci)(2T+ 1)^ ^ Mjv, + (^) = c 1 exp{Nc 1 ){2T + 1) N Z Nt+ . □ 


Lemma 10.4. There exists C 3 > 0 such that 

Zn- 1 ,++ 


■'N 


< C3 exp(iVc3) • N 


-26 N 


N = 1 , 2 ,. 


Proof. We have 
Zn 


(146) N 


Zn- 1 , 


++ 


W E 


M»(() 

Zn- 1 , 


N 

EE 

k= 1 


i#E++(* (fc) ) 


++ 


4 fe) -e N—i 


e n r (| m ~1! l +1) r ( ' m :. < n +(,) exp (- iv 1 /(^ 


o(P) 1 


Zn—i, 


++ 


„ = ,<» a r(|m - 4‘>|)r(|m - if >I + 1-9) 


N 


where £^ varies over ( |139[ ) with N replaced by N — 1, and we use the notation = — 00 , 

/ iy\ _ Q 

t N = + 00 . Note that by the definitions, the sum k m „ is always non-empty. Take £ 


m =4-i+ 0 


satisfying (136) with N replaced by N — 1 and let s^ k \£) denote the corresponding £' of Lemma 


10.2 Then (146) implies 


Zn 


n 


(147) AT — >(2T+l)- jV ^^ 

Zn- 1 ,++ 


j( fc ) 

%-!,++ 


(»“>(<)) 


k= 1 


Zn~ 


N- 1,++ 


9 E 1 r(|m - s (k \£)j\ + l)r(|m - s^(£)j\ + 0 ) / nn 

f-i r d m -s( k \£)i\)T(\m-s( k )(£)i\ + 1-6) V 

m=s^>(Z) h _ 1 +6 *—r 


N 


where £ varies over (136). 
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Using Stirling’s formula, its corollary (91), and definition of the state space (139), we see that 
the product of Gamma functions in the above formula is bounded from below by 

exp(— Ng) _ N 2 e N 
c 4 


with a constant C 4 > 0. Let v be the maximum of V(x) over [—2T, 2T], Then (147) and Lemma 

n sW(t%-e 


10.2 imply that for C 5 > 0 


N Z N > N 26 N exp(-c 5 lV) 


Zn- 


1 ,++ 


C5 


I-2T<m<2T exp ( — Nv) 

k =1 t m=s( k )(t) k _ 1 +e 


Note that for every i satisfying (136) there exists at least one k, such that for at least one 
m G s( k \£)k-i + 6, s^ k \i)k -1 + 0 + 1 ,..., s( k \l)k — 9, we have \m\ < 2 T (here again we use the 
notation = — 00 , s^(£)n = + 00 ). Therefore, 

n J n > ++ (f) expf-JVp) = N™« eXp( - C5JV) exp(-iVp). □ 

Zn- 1,++ c 5 i ^ 

Lemma 10.5. There exists C 4 > 0 such that 

Zn. 


c 5 


(148) 


< C 4 exp(c4^7), IV = 1,2, 


z n,++ 

Proof. Dehne the random probability measure vn through 

1 N 

UN = Jj ^2 S ti/(N+I)i {fit ■ ■ -An) is Pjv,++-distributed 


i =1 


Then 

(149) 


Zn,- 

Zn, 


++ 


E 


MivMg 

Zn. 


= Ep 


++ 


JV,++ 


exp ^ N J V(x)vN(dx] 


In order to bound (149) we start with a lower bound 


Z N ,++ > M N)++ (0, 26,..., N6) > exp(-c 5 fV 2 ) • N [ 


2617V 2 


for some C 5 > 0. On the other hand, using (91) we obtain 
lnMjv, ++ (£i,... ,£n) 


< exp 26 N 2 In N - N 2 




x^y 


V{x) + V(y) 


- 6ln\x - y\ ) VN(dx)v N {dy) + cqN 2 


for some cq > 0. Assumption (134) and inequality In \x — y\ < In(|rc| + 1) + ln(|y| + 1) imply that 
that there exist C 7 > 0 such that for all x / y, 


V(x) + V(y) — 26 In \x — y\ 

> ^- £ (V{x) + v (y)) + 


V{x) 

1 -\- S 


- 201 n(|a;| + 1 ) ) + 


( v(y) 

\ 1 + £ 


> 


— 29 ln(|y| + 1 ) 


(U(x) + V{y)) - c 7 . 


1 + £ 
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Therefore, for each L > 0 we have 


Ep 


JV.++ 


V(x)/j,N(dx) > L 


< exp ( c 5 N 2 + c e N 2 + ^-N 2 - ^ 


N 


w 2 i)n e “p 

2—1 \mEZ+20 


2(1 + £) 


NV (m/N) 


Hence, there exists cs > 0, and Lq such that for all L > Lq, 


Ep 


N,++ 


V(x)fiN(dx) > L 


< exp (—c%LN 2 ) 


Together with (149) this implies (148). 


□ 


Proof of Theorem 10.1. We have 

Z N- 1,+ 


(150) W N (l N > DN) = 


Z N 




N -1 

n 

2—1 _ 


m> max(£jv_i, DN) 

T(m — Ij + l)T(m — lj + 0) 
T(m - Zi)T(m - £ + 1 - 0) 


ex P MM ^ 


where t ranges over (|136|) with N replaced by IV — 1. The combination of Lemmas 

< C 5 exp(lVc 5 ) • N~ 29N . 


implies that for C 5 > 0 we have a bound 

Zn-i- 


10.3 


and 


10.5 


Z N 


For the product of gamma-functions in (150) we use Stirling’s formula and its corollary (91), 
which yields 


N—l 

n 

2—1 


r(m -£j + 1 )T(m -ij + 6) 

T(m - ii)T(m -ii + l-0) 


N-l 


m. ii 
N~ N 


26 


< c 6 exp(c 6 N)N 26N J] 

2 — 1 

< c 6 exp(c 6 iV)iV 20 JV 2 (20+1)(iV - 1) Yl 


N -1 


i =1 . 


26 


+ 


u 

N 


26 


It follows that for c 7 > 0, 

(151) P n (£ n > DN) < c 7 exp(c 7 lV)^P 7 v- 1 ,+(^) exp ( -NV 


N-l 

n 

2—1 _ 


26 


+ 


N 


26 > 


exp u 


Further, take <5 > (20) , and observe that when D > H + 1, 


£ exp 

m>max(^iv— i,DN) 


00 00 

(-5H(^))< £ exp(-2«J01n(m/JV)) = N 2Se £ 


m=|ZWJ 
< 2ND 1 ~ 269 . 


m=[DN} 
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Also due to (134), when m > DN > (H + 1 )N, we have 


29 


+ 


N 


29 > 


exp 


Therefore, writing N(V(^)) in (151) as 

ATIVI^!I = 1 + £/2 



< const 


N 

m 


Qe 


< const • D 


—Oe 






and noticing that the last term is positive when m > DN and N is large, we conclude that for 
some c 8 > 0, all D > H + 1 and all N > Nq 

P n (£ n > DN) < c 8 exp(c 8 IV) • N 2 • D~ 9eN . 

Choosing D large enough, we obtain the desired exponential estimate for P/v(C/v > DN). The 
estimate for P/v(^i < —DN) is obtained in the same way — the only difference is that we now 


need to bound 


7 ( 1 ) 

Z N 


instead of 


'^N- 1 ,+ 

Z N ■ 


but for that we use Lemma 


10.4 


□ 


Remark 10.6. It is very plausible that one can similarly establish an analogue of Theorem 10.1 


for more general models in the framework of Section [ 3 ] with a±(N) = —00 and b}.(N) = + 00 . 
The only necessary modification in the above proofs is in Lemma 10. 2[ where we should take into 
account that the Lipshitz property of V(x) might fail near the endpoints a/(IV), bi(N), cf. (72). 


We will not address here the exact conditions on V(x) under which an analogue of Lemma 10.2 
holds for the models in framework of Section [H 
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